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The cubic complex Ginzburg-Landau equation is one of the most-studied nonlinear equations in 
the physics community. It describes a vast variety of phenomena from nonlinear waves to second- 
order phase transitions, from superconductivity, superfluidity and Bose-Einstein condensation to 
liquid crystals and strings in field theory. Our goal is to give an overview of various phenomena 
described the complex Ginzburg-Landau equation in one, two and three dimensions from the point 
of view of condensed matter physicists. Our approach is to study the relevant solutions to get an 
insight into nonequilibrium phenomena in spatially extended systems. 
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I. PRELIMINARY REMARKS 
A. The equation 

The cubic complex Ginzburg-Landau equation (CGLe) is one of the most-studied nonlinear equations in the physics 
community. It describes on a qualitative, and often even on a quantitative level a vast variety of phenomena from 
nonlinear waves to second-order phase transitions, from superconductivity, superfluidity and Bosc-Einstoin conden- 
sation to liquid crystals and strings in field theory (Kuramoto, 1984, Cross and Hohenberg, 1993, Newell et al, 1993, 
Pismen, 1999, Bohr et al, 1998, Dangelmayr and Kramer, 1998). 

Our goal is to give an overview of various phenomena described by the CGLe from the point of view of condensed 
matter physicists. Our approach is to study the relevant solutions to get insight into nonequilibrium phenomena in 
spatially extended systems. More elementary and detailed introductions into the concepts underlying the equation 
can be found in Manneville (1990), van Saarloos (1993), van Hecke et al (1994), Nicolis (1995) and Walgraaf (1997). 

The equation is given by 

dtA = A + {l + ib)AA-{l + ic)\A\''A, (1) 

where A is a complex function of (scaled) time t and space x (often in reduced dimension I? = 1 or 2) and the 
real parameters h and c characterize linear and nonlinear dispersion. The equation arises in physics in particular as 
a "modulational" (or "envelope" or "amplitude" ) equation. It provides a reduced, universal description of "weakly 
nonlinear" spatio-temporal phenomena in extended (in x) continuous media whose linear dispersion is of a very 
general type (sec below) and which are invariant under a global change of gauge (multiplication of A by exp(i$)). 
This symmetry typically arises when A is the (slowly varying) amplitude of a phenomenon that is periodic in at least 
one variable (space and/or time) as a consequence of translational invariancc of the system. 

The assumptions of slow variation and weak nonlinearity are valid in particular near the instability of a homogeneous 
(in x) basic state and Eq. (1) can be viewed as a (generalized) normal form of the resulting "primary" bifurcation. 
Then, in analogy with phase transitions, A is often called an order parameter. 
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To see more clearly the analogy with the order parameter concept we write the equation in the unsealed form, ^ as 
derived for example from the underlying set of basic (e.g. hydrodynamic) equations for a definite physical situation 



T{dfA - Vg ■ Vi) = e(l + ia)A + ^^(1+ ib)AA - g{l + ic)\A\ 



'A, 



(2) 



Equation (1) is obtained from (2) by the transformations A = {e/gY^^A exp{—i{ea/T)f), i = {T/t)t, and x — Vgi = 
(^/e^/^)x. The case e > and g > was assumed. Otherwise the signs in front of the first and/or last term on the 
right-hand side of (1) have to be reversed. The physical quantities u(i, f) (temperature, velocities, densities, electric 
field etc.) are given in the form 



(c.c.= complex conjugate, h.o.t.= higher-order terms). If the phenomena occur in (thin) layers, on surfaces, or 
in (narrow) channels, then U;, derived from the linear problem, describes the spatial dependence of the physical 
quantities in the transverse z direction(s). U; is a linear eigenvector and Wc'fc the corresponding eigenvalues. In the 
case of periodically driven systems, U; would include a periodic time dependence. 

In order to identify the character of the various terms in the linear part of (2) one may also consider the dispersion 
relation obtained from (2) and (3) for small harmonic perturbations of the basic state A = 



Here q is the wavevector in the physical system and \ = <j — ito\s the complex growth rate of the perturbation, r is a 
characteristic time, ^ is the coherence length, Vg a linear group velocity, and ea/r a correction to the Hopf frequency 
ujc- e measures in a dimensionless scale the distance from threshold of the instability, i.e. e = {R — Rc)/ Rc, with R 
the control parameter that carries the system through the threshold at Rc- Note that there is an arbitrary overall 
factor in Eqs. (2) and (4) which is fixed by the ultimately arbitrary choice of the definition of e. The value of the 
nonlinear coefficient g in Eq. (2) depends on the choice of the normalization of the linear eigenvector U;. 

Now we can proceed to summarize the conditions for validity of the CGLe. The following four points are to some 
extent interrelated. 

a) Correct choice of order parameter space, i.e. a single complex scalar: First of all this necessitates that Uc 
and/or qc are nonzero, because otherwise one would expect a real order parameter as in simple phase transitions. An 
exception is the transition to superconductivity and superfluidity where the order parameter is complex for quantum 
mechanical reasons (see below). Moreover, if qc = 0, one may run into problems with conservation laws which 
frequently exclude a homogeneous change of the system. In this case of long-wavelength instabilities often somewhat 
different order-parameter equations arise (see e.g. Nepomnyashchii, 1995a). Secondly, a discrete degeneracy (or near 
degeneracy) of neutral modes is excluded, which may arise by symmetry (see below for an example) or by accident. 
If the eigenvectors of the different modes are different one would need several order parameters and a set of coupled 
equations. If the eigenvectors coincide (or nearly coincide), which may happen at (or near) a co-dimension-2 point, 
one can again use one equation, which would now contain higher space or time derivatives. 

b) Validity of the dispersion relation (4): since there is no real contribution linear in q — qc and since is a 
positive quantity the real growth rate a has a minimum at qc- In more than ID this excludes an important class of 
systems, namely isotropic ones with gc 7^ like Rayleigh-Benard convection in simple 2D fluid layers. There one has 
a continuous degeneracy of neutral linear modes. The neglect of terms of higher order in e and in ^— is usually 
justified near the bifurcation. 

c) Symmetries: translation invariance in x and t. Actually the CGLe incorporates translational invariance with 
respect to space and/or time on two levels. One is expressed by the global gauge invariance, which in the CGLe can 
be absorbed in a shift of x and/or time t. Note that this invariance excludes terms that are quadratic in A. The other 
is expressed by the autonomy of the CGLe (no explicit dependence on space and time) . The two invariances reflect 
the fact that the fast and the slow space and time scales are not coupled in this description. This is an approximation 
which cannot be overcome by going to higher order in the expansion in terms of amplitude and gradients. The coupling 
effects are in fact nonanalytic in e ( "non-adiabatic effects", see, e.g. Pomeau (1984), Kramer and Zimmermann (1985), 
Bensimon et al (1988)). 



^except maybe for a simple rescaling and rotation of the coordinate system, see below 



u = ie'(«<=-^-'^=*'U;(5} + ex. + h.o.t. 



(3) 



tA = -iTLJc + ivg ■ {q- qc) + e(l -|- ia) - ^^(1 -h ib){q- g^)^. 



(4) 
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d) Validity of the (lowest-order) weakly nonlinear approximation: we will deal mostly with the case of a supercritical 
("forward", or "normal") bifurcation where g > and then higher-order nonlinearities in Eq. (1) can be neglected 
sufficiently near threshold. If the nonlinear term in Eq. (1) has the opposite sign, which corresponds to a subcritical 
("backward" or "inverse") bifurcation, higher-order nonlinear terms are usually essential. However, even in this case, 
there exist for sufficiently large values of |c| relevant solutions that bifurcate supercritically, which will be discussed 
in Sec. VI.A.l. 

From the linear theory we can now distinguish three classes of primary bifurcations where Eq. (1) arises: 

i) u>c = 0, gc 7^ 0: for such stationary periodic instabilities A is real, and in fact all the imaginary coefficients 
(including the group velocity Vg) vanish. Generically, reflection symmetry is needed (see below). Equation (1) then 
reduces to the "real" Ginzburg-Landau Equation (GLe) 

dtA = A + AA- \AfA (5) 

which one might also call the "Complex Nonlinear Diffiision Equation" in some analogy with the Nonlinear 

Schroedinger Equation (see below). Examples that display such an instability are Rayleigh-Bcnard convection in 
simple and complex fluids, Taylor-Couette flow, electroconvection in liquid crystals and many others. In more than 
ID there is the restriction mentioned under b). Thus in isotropic 2D systems the dispersion relation is changed and 
the Laplacian in Eq.(5) has to be substituted by a different differential operator. The corresponding equation derived 
by Newell and Whitehead (1969) and by Segel (1969) was in fact the first amplitude equation that included spatial 
degrees of freedom. It is applicable only for situations with nearly parallel rolls, which is in isotropic systems an 
important restriction. 

So in more than one dimension the system must be anisotropic, which is the case in particular for convective 
instabilities in liquid crystals (Kramer and Pesch, 1995), but holds also for Rayleigh-Benard convection in an inclined 
layer (Daniels et al, 2000) or in a conducting fluid in the presence of a magnetic field with an axial component (Eltaycb, 
1971). Also the Taylor-Couette instability in the small gap limit can be viewed as an anisotropic quasi-2D system. In 
2D Eq. (5) was first considered in the context of electrohydrodynamic convection in a planarly aligned nematic liquid 
crystal layer (Pesch and Kramer, 1986, Bodenschatz et al, 1988a, Kramer and Pesch, 1995, for review see also Buka 
and Kramer, 1996) The Laplacian in Eq. (2) is obtained after a linear coordinate transformation. 

ii) ujc 7^ 0, ^ = 0: The prime example for such oscillatory uniform instabilities are oscillatory chemical reactions 
(see e.g. de Wit, 1999). In lasers (or passive nonlinear optical systems) this type may also arise (Newell and Moloney, 
1992). In hydrodynamic systems such instabilities are often suppressed by mass conservation (see, however, Borzsonyi 
et al, 2000). Isotropy does not cause any problems here and the Laplacian applies directly. In the presence of 
reflection symmetry the group velocity term in Eq. (2) is absent. The spatial patterns obtained from A reflect 
directly those of the physical system. The imaginary parts proportional to b and c pertain to linear and nonlinear 
frequency change (renormalization) of the oscillations, respectively. In most systems the nonlinear frequency change 
is negative (frequency decreases with amplitude), so that c < with our choice of signs. Coefficients for the CGLe 
have been determined e.g. from experiments on the Belousov-Zhabotinsky (BZ) reaction (Hynne et al, 1993, Kramer 
et al, 1994). 

iii) uJc ^ Q,qc ^ 0: This oscillatory periodic instability occurs in hydrodynamic and optical systems. The best- 
studied example is Rayleigh-Benard convection in binary mixtures, although here the bifurcation is in the accessible 
parameter range mostly subcritical (Schopf and Zimmermann, 1990, Liicke et al, 1992). Also, in 2D, the system 
is isotropic, so that the simple CGLe (2) is not applicable. Other ID examples are the oscillatory instability in 
Rayleigh-Benard convection in low Prandtl number fluids, which in 2D occurs as a secondary instability of stationary 
rolls. In a ID geometry with just one longitudinal roll it can be treated as a primary bifurcation (Janiaud et al, 1992). 
Other examples include the wall instability in rotating Rayleigh-Benard convection (Tu and Cross, 1992, van Hecke 
and van Saarloos, 1997, Yuanming and Ecke, 1999) and hydrothermal waves, where the coefficients of the CGLe were 
determined from experiment (Burguette et al, 1999). In 2D the prime example is the electrohydrodynamic instability 
in nematic liquid crystals in thin and clean cells (otherwise one has the more common stationary rolls) (Trciber and 
Kramer, 1998). Actually in such an anisotropic 2D system one is lead to a generalization of Eq. (1) where the term 
ibAA is replaced by a more general bilinear form iipid^ + b2dy), see Sec. VI. D. 

Most of the oscillatory periodic systems just mentioned have reflection symmetry and then one has to allow for 
the possibility of counter-propagating waves which makes a description in terms of two coupled CGLes necessary (see 
Cross and Hohenberg, 1993). The degeneracy between left and right traveling waves can be lifted by breaking the 
reflection symmetry by applying additional flelds or an additional flow. In this situation one roll system is favored over 
the other and, if the effect is strong enough, a single CGLe can be used. Breaking reflection symmetry in stationary 
periodic instabilities (case i)) the rolls will generically start to travel and one indeed arrives at an oscillatory periodic 
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instability (case iii)). This has been studied experimentaUy by applying a through flow in thermal convection (Poehcau 
and Croquette, 1984) or in the Taylor-Couette system (Tsameret and Steinberg, 1994; Babcock, Ahlers, and Cannell, 
1991) or by non-symmetric surface alignment in electroconvection of nematics (pretilt or hybrid alignment, see e.g. 
Krekhov and Kramer, 1996). 

Since the drift introduces a frequency, for sufficiently strongly broken reflection symmetry the distinction between 
cases i) and iii) is lost, as is obvious in open-flow systems (Leweke and Provansal, 1994, 1995; Roussopoulos and 
Monkewitz, 1996). 

The CGLe may also be viewed as a dissipative extension of the conservative nonlinear Schrodinger equation (NLSe) 

idtA = AA± \A\^A , (6) 

which describes weakly nonlinear wave phenomena (Newell, 1974). The prime examples arc waves on deep water, 
(Dias and Kharif, 1999) and nonlinear optics (Newell and Moloney, 1992). The conservative limit of Eq. (1) is 
obtained by letting in Eq. (2) Vg, e, ^, 5 — > with ^^6, gc remaining nonzero, so that b, c — > oo, and rescaling space 
and the amplitude. 

B. Historical remarl<s 

Four key concepts come together in the CGLe philosophy: 

• Weak nonlinearity, which amounts to an expansion in terms of the order parameter \A\. This concept goes 

back to Landau's theory of second-order phase transitions (Landau, 1937a). Landau also employed this type of 
expansion in his attempt to explain the transition to turbulence (Landau, 1944). In the context of stationary, 
pattern-forming, hydrodynamic instabilities the weakly nonlinear expansion leading to a solvability condition at 
third order was introduced by Gorkov (1957) and Malkus and Veronis (1958). One should also mention the work 
of Abrikosov in 1957 (for review see Abrikosov (1988)) where he presents the theory of the mixed state of type 
II superconductors in a magnetic field based on the Ginzburg-Landau theory of superconductivity. The mixed 
state is a periodic array of flux lines (or vortices) corresponding to topological defects (see below). Abrikosov 
introduced a weakly nonlinear expansion to describe this state valid near the upper critical field. 

• Slow relaxative time dependence was first used by Landau in the above-mentioned paper on turbulence in 1944. 
In the context of pattern-forming instabilities it goes back to Stuart (1960). 

• Slow nonrelaxative time dependence with nonlinear frequency renormalization in the complex-amplitude formu- 
lation was introduced by Stuart in 1960 using multi-scale analysis. Of course perturbation theory for periodic 
orbits (in particular conservative Hamiltonian systems) is a classical subject that was treated by Bogoliubov, 
Krylov and Mitropolskii in 1937 (for review, see Bogoliubov and Mitropolskii, 1961). 

• Slow spatial dependence was included already by Landau (1937b) in the context of X-ray scattering by crystals 
in the neighborhood of the Curie point. However, the concept became known with the success of the (stationary) 
phenomenological Ginzburg-Landau (GL) theory of superconductivity [Ginzburg and Landau (1950)]. 

The stationary GL theory for superconductivity has a particular resemblance to the modulational theories of pattern- 
forming systems because the order parameter is complex, although for a very different reason. Superconductivity being 
a macroscopic quantum state requires an order parameter that has the symmetries of a wave function. In spite of the 
different origin one has many analogies. However, in superconductors the time dependence is rendered nonvariational 
primarily through the coupling to the electric field due to local gauge invariance (see, e.g. Abrikosov, 1988), a 
mechanism that has no analog in pattern- forming systems. 

The time-dependent GL theory for superconductors was presented (phenomenologically) only in 1968 by Schmid, 
(derived from microscopic theory shortly afterwards by Gorkov and Eliashberg (1968)), when the first modulational 
theory was derived in the context of Rayleigh-Benard convection by Newell and Whitehead (1969) and Segel (1969). 
Eq. (5) with additional noise term has been studied intensively as a model of phase transitions in equilibrium systems, 
see e.g. Hohcnberg and Halperin (1977). 

The full CGLe was introduced phenomenologically by Newell and Whitehead (1971). It was derived by Stewartson 
and Stuart (1971) and DiPrima, Eckhaus and Segel (1971) in the context of the destabilization of plane shear flow, 
where its applicability is limited by the fact that one deals with a strongly subcritical bifurcation. In the context of 
chemical systems the CGLe was introduced by Kuramoto and Tsuzuki (1974). 
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There exists an extended mathematical hteraturc on the CGLe, which we will touch rather little, see e.g. Doering 
et al (1987, 1988), Levermore and Stark (1997), Doelman (1995), van Harten (1991), Milke and Schneider (1996), 
Schneider (1994), Melbourne (1998), Milke (1998). 

C. Simple model - vast variety of effects 

Clearly the CGLe (1) may be viewed as a very general normal- form type equation for a large class of bifurcations 
and nonlinear wave phenomena in spatially extended systems, so a detailed investigation of its properties is well 
justified. The equation "interpolates" between the two opposing limits of the conservative NLSe and the purely 
relaxative GLe. The CGLe world lies between these limits where new phenomena and scenarios arise, like sink and 
source solutions (spirals in 2D and filaments in 3D), various core and wave instabilities, nonlinear convective versus 
absolute instability, screening of interaction and competition between sources, various types of spatio-temporal chaos 
and glassy states. 

II. GENERAL CONSIDERATIONS 

In this Section we will study general properties of Eq. (1) relevant in all dimensions. 

A. Variational case 

For the case of 6 = c it is useful to transform into a "rotating" frame A — * Aexp(i6f). Then Eq. (1) goes over into 

dtA ={1 + ih) {A + I^A- \A\'^A) . (7) 
Eq. (7) can be obtained by variation of the functional 

V = jud''r, U = -\Af + ^\A\''+\VAf. (8) 

leading to dtA = -(1 -|- ib)6V/5A* and 

^^^=-YT^ J IdtAl'd'^r. (9) 

One sees that for all non-infinite b the value of V decreases, so the functional (8) plays the role of a global Lyapunov 
functional or generalized free energy (V is bounded from below). The system then relaxes towards local minima of 
the functional. In particular the stationary solutions A(f) of the GLe (b=0) correspond in the more general case to 
j\^g-tbt corresponding stability properties. 

In the NLSe limit 6 ^ oo the functional becomes a Hamiltonian, which is conserved. More generally, the NLSe is 
obtained from the CGLe by taking the limit 6, c ^ oo without further restrictions. After rescaling one obtains Eq. 
(6). One sees that the equation comes in two variants, the focusing (+ sign) and defocusing (—sign) case (the notation 
comes from nonlinear optics. In ID it is completely integrable (Zakharov and Shabat, 1971). In the focusing case 
it has a two-parameter family of "bright" solitons (irrespective of space translations and gauge transformation). In 
D> 1 solutions typically exhibit finite-time singularities (collapse) (Zakharov, 1984, for recent review see Robinson, 
1997). In the defocusing case one has in ID a three-parameter family of "dark solitons" that connect asymptotically 
to plane waves and vortices in 2 and 3D. We will here not discuss the equation since there exists a vast literature on 
it (see e.g. Proceedings of Conference on The Nonlinear Schrodinger Equation,1994). It is useful to treat the CGLe 
in the limit of large b and c from the point of view of a perturbed NLSe. 

B. The amplitude-phase representation 

Often it is useful to represent the complex function A by its real amplitude and phase in the form A = Rexp{i6). 
Then Eq. (1) becomes 
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dtR = [A - {Vef]R - h{2Ve ■ Vi? + RAO) + (1 - R^)R , (10) 
Rdte = b[A - {Vef]R + 2V6» • VR + RAO - cR^ . 

For 6 = this corresponds to a class of reaction-diffusion equations called A — w systems, which are generally of the 
form 

dtR-^[A~ {Ve)^]R + RX{R) , (11) 
Rdt0 = 2Ve ■ VR + RAO + Ru{R) . 

Such equations have been studied in the past extensively by applied mathematicians (see Hagan, 1982, and references 
therein). Clearly, one can combine Eqs. (10) in such a way that the right-hand sides are those of a A — w system. 



C. Transformations, coherent structures, similarity 

The obvious symmetries of the CGLe are time and space translations, spatial reflections and rotation, and global 
gauge (or phase) symmetry A Ae^'^. The transformation A,h,c^ A* ,—h, — c leaves the equation invariant so that 
only a half plane within the &, c parameter space has to be considered. 

Other transformations hold only for particular classes of solutions. To see this it is useful to consider the following 
transformation 

A{x,t) = e'W-^-'^*)B(f - vt,t), (12) 

leading to 

5tB = [ct + w5 • V + (1 + ih)A - (1 + ic)\B\^]B (13) 

with (T = 1 + icj — (1 + ib)Q'^ , w = v — i{l + ih)Q. Most known solutions are either of the coherent- strucAure type, 
where B depends only on its first argument (i.e. with properly chosen uj it is time independent in a moving frame), 
or are disordered in the sense of spatio-temporal chaos. Coherent structures can be localized or extended. The "outer 
wavevcctor" Q could be absorbed in B (then u) would be replaced by cj — Qv) . It may be useful to introduce Q when 
the gradient of the phase of B, integrated over the system, is zero (or at least small). 

With a little bit of algebra one can now derive a useful similarity transformation that connects coherent structures 
along the lines (6 — c)/(l + he) = const in parameter space. Defining 

f=Pr', |B|=7|B'|, (14) 

the transformation relations between unprimed and primed quantities can be written as 

b-c _ b'-c' 
1 + bc ~ 1 + b'c' ' 

The relations (15) are independent of the nonlinear part of (1) and therefore survive generalizations. Solutions 
with V = remain stationary (for b < oo) and that transformation was given by Hagan, 1982. Clearly one can very 
generally transform to 6 = 0, where the CGLe represents a \ — uj system. Note that for v one cannot have Q = 
and Q' = 0. The similarity line b = c (vanishing group velocity, see below) includes the real case. 

Note, that the stability limits of coherent states are in general not expected to conform with the similarity transform. 
By taking in Eqs. (16) the limit b oo, one finds c' = c. Then the 1 in the factor 1 + ib can be dropped and changes 
in b can be absorbed in a rescaling of length. The similarity transformation then connects solutions with arbitrary 
velocity, which is a manifestation of a type of Galilean invariance (van Saarloos and Hohenberg, 1992). Thus in this 
limit solutions appear as continuous families moving at arbitrary velocity. 

Similarly, by letting c ^ oo, 6 tends to a constant. Then, in addition to the 1 in the factor 1 +ic, the linear growth 
term can be dropped, and then changes in c can be absorbed in a rescaling of B as well as length and time. The 



1 + bc 1 + c'^ 
1 + b'c' 1 + c2 



(16) 
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similarity transformation then turns into a scaling transformation. Thus in this limit solutions appear as continuous 
families of rescaled functions. 

For b and c — > oo one has both transformations together, so solutions appear generically as two-parameter families. 
Indeed, one is then left with the NLSe. 

D. Plane-wave solutions and their stability 

The simplest coherent structures are the plane-wave solutions 

A=^/l^^exp[i{Q■r-u)p{Q)t + ^)], F^ = l-Q^, u}p{Q) = c{l - Q^) + bQ^ . (17) 

{(j) is an arbitrary constant phase) which exist for < 1. To test their stability one considers the complex growth 
rate A of the modulational modes. One seeks the perturbed solution in the form 

A = {F + Sa+ exp[Ai -h ik • r] -|- Sa- exp[A*i - ik • r]) exp[i(Q • r - tot)] (18) 

where k is a modulation wave vector and 6a± are the amplitudes of the small perturbations. One easily finds the 
expression for the growth rate A (Stuart and DiPrima, 1980): 

+ 2{F'^ + 2ibQ-\<i + k^)X 

+(1 + 62)(fc* - 4(Q • k)4) + 2F^[{1 + bc)k^ + 2i{b - c)Q • k] = 0. (19) 
By expanding this equation for small k one finds 

A = -iVgk - £>2fc^ + i^gk^ - D^k'^ + 0{k^), (20) 

with 

Vg = 2(6 - c)Qk 
D. = l + bc-'^^±^ 



1-Q2 

2[b{l-Q^)-2cQl]{l+d')Qk 

(1-02)2 



= 2(1 -qV [^'(1 - Q'f - 12&c(l - Q')Ql + 4(1 + 5c')Qt] 



(21) 



where = Q • k is the component of Q parallel to k. The quantities Vg, D2, i^g and D4 for k || Q will be denoted by 

V^ll , Z?2|| , ^g\\ and D411 . Similarly, for k _L Q we may use the notation Vgj_, D2_l, ^gi_ and D^^. Clearly the longitudinal 
perturbations with k || Q are the most dangerous ones. The solutions (17) are long-wave stable as long as the phase 
diffusion constant i32|| is positive. Thus one has a stable range of wave vectors with < Q\ = (1 + 6c)/(3 + 2(? + he) 
enclosing the homogeneous [Q = 0) state as long as the "Bcnjamin-Feir- Newell criterion" 1 + 6c > holds. This 
criterion conforms with the similarity transform (16). The condition D2\\ > is the (generalized) Eckhaus criterion. 
For 6 = c it reduces to the classical Eckhaus criterion < Q\ = 1/3 for stationary bifurcations. We will call 
the qTiadrants in the 6, c plane with 6c > the "defocusing quadrants" . Otherwise we will speak of the "focusing 
quadrants" . 

From Eqs. (19,20) one sees that for 6 — c ^ and Q 7^ the destabiHzing modes have a group velocity Vg = V qVL = 
2(6 — c)Q, so the Eckhaus instability then is of convective nature and does not necessarily lead to destabilization of 
the pattern (see next subsection). 

The Eckhaus instability signalizes bifurcations to quasiperiodic solutions (including the solitary limit), which are 
of the form (13) with periodic function B{x — Vt). It becomes supercritical before the Benjamin-Feir-Ncwell criterion 
is reached (Janiaud et al, 1992) and remains so in the unstable range. The bifurcation is captured most easily in the 
long- wave limit by phase equations, see below. A general analysis of the bifurcating solutions has recently been done 
(Brusch et al 2000), see subsection III.C. 

It is well known that the Eckhaus instability in the CGLe is not in all cases of long-wave type. Clearly, a necessary 
condition for it to be the case is that Dm is positive where D2\\ changes sign. 
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From the above expressions one deduces that this fails to be the case for \b\ > 64(c), where 
64 = 



2|c|(l - c^) - V4c2(l - c2)2 - (1 - 3c2)(l + 5c2) /(l-3c^), 600. (22) 

Thus the range is in the defocusing quadrant, far away from the Benjamin- Feir-Newell stabihty hmit. We are presently 
not aware of effects where this phenomenon is of relevance. 

E. Absolute versus convective instability of plane waves 

For a nonzero group velocity Vg — Vgfi — 2(6 — c)Q the Eckhaus criterion can be taken only as a test for convective 
instability. In this case a localized ID initial perturbation 6*0(2;) of the asymptotic plane wave, although amplified 
in time, drifts away and does not necessarily amplify at a fixed position (Landau and Lifshitz, 1959). For absolute 
instability localized perturbations have to amplify at fixed position. The time evolution of a localized perturbation is 
in the linear range given by 

/oo 
dfc/(27r)S'o(fc) exp{ikx + X{k)t) (23) 
-00 

where S'o(fc) is the Fourier transform of So{x). ^ The integral can be deformed into the complex k-plane. In the limit 

t ^ 00 the integral is dominated by the largest saddle point ko of \{k) (steepest descent method, see e.g., Morse and 
Feshbach, 1953) and the test for absolute instability is 

i?e[A(fco)] > with dkX{ko) = (24) 

The long-wavelength expansion (20) indicates that at the Eckhaus instability, where D2 becomes negative, the system 
remains stable in the above sense. When D2 vanishes and Q the main contribution comes from the term linear 
in k that then can suppress instability. 

In the following, results of the analysis arc shown in the b—c plane (Aranson et al, 1992; Weber et al, 1992). In Fig. 1 
the scenario is demonstrated for four cuts in the b, c, Q space. The stable region (light gray) is limited by the Eckhaus 
curve, which terminates at Q = on the curve. To the right of it there exists a convectively unstable wavenumber 
band < Qai < \Q\ < Qa2 (dark grey). Qai goes to zero on the Benjamin- Feir-Newell curve as (1 -|- bc)^/^, which 
can be seen from the long-wavelength expansion (20) with the fourth-order term included. Moving away from the 
Benjamin- Feir-Newell line (into the unstable regime) Qai increases and Qa2 decreases until they come together in a 
saddle-node-type process at a value Qac- Beyond Qac there are no convectively unstable plane waves. The saddle 
nodes Qac are shown in the b, c plane in Fig. 2 (curve Al). Thus, convectively unstable waves exist up to this curve. 
The Benjamin- Feir-Newell criterion curve up to which convectively unstable waves exist is also included. The other 
curves include the Eckhaus instability and the absolute stability limit for waves with wavenumber selected by the 
stationary hole solutions (see Sec. III.B). 

Eq. (24) by itself gives only a necessary condition for absolute stability. However, it is also sufficient, as long as 
one of the two roots fci,2 of the dispersion relation, which collide at the saddle point when Re[X] is decreased from 
positive values to zero, say fci, is the root that produced the convective instability, and if ^2 does not cross the real 
A;— axis before ki does, when Re[X] is increased from zero. (Note that ki has to cross the axis when Re[X] is increased. 
Note also that, because of the symmetry X,k ^ A*, — fc* one has parallel processes with opposite sign of /to[A]). For 
a discussion of the underlying "pinching condition", see e.g. Brevdo and Bridges (1996). In the parameter range 
considered here the sufficient condition is fulfilled. 

Clearly the absolute stability boundary can be reached only if reflection at the boundaries of the system is sufficiently 
weak, so in general one should expect the system to lose stability before the exact limit is reached. In many realistic 
situations in the CGLe the interaction of the emitted waves with the boundaries leads to sinks (shocks) in analogy 
to the situation where different waves collide. These shocks are strong perturbations of the plane wave solutions but 
they absorb the incoming perturbations. Thus here, even for periodic boundary conditions, the absolute stability 
limit is relevant. Moreover, in infinite systems, states with a cellular structure made up of sources surrounded by 



^It can be shown strictly that the destabilization occurs at first for purely longitudinal perturbations Q || k 
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sinks can exist in the convectively unstable regime. However, the convectively unstable states are very susceptible to 
noise, which is exponentially amplified in space. The amplification rate goes to zero at the convective stability limit 
and diverges at the absolute stability limit. 

The concept of absolute stability is relevant in particular for the waves emitted by sources and possibly also for some 
characteristics of spatio-temporally chaotic states. When boundaries are considered one also has to allow for a linear 
group velocity term in the CGLe (see Eq. (2)). Such a term does not change the convective stability threshold, but 
clearly the absolute stability limit is altered, and this is important in particular in the context of open flow systems. 
The saddle point condition in (24) ensures the existence of bounded solutions of the linear problem that satisfy 
nonperiodic boundary conditions (irrespective of their precise form) at well separated side walls (sometimes called 
"global modes"), because for that purpose one needs to superpose neighboring (extended) eigenmodes, which are 
available precisely at the saddle point. This is an alternative view of the absolute instability (Huerre and Monkowitz, 
1990, Tobias and Knobloch, 1998, Tobias, Proctor, and Knobloch, 1998). 

Often the condition (24) coincides with the condition that a front invades the unstable state in the upstream 
direction according to the linear front selection criterion (marginal stability condition, see, e.g. van Saarloos (1988)). 

F. Collisions of plane waves and effect of localized disturbances 

The nonlinear waves discussed above have very different properties from linear waves. In particular, when two waves 
collide they almost do not interpenetrate. Instead a "shock" (sink) is formed along a point (ID), line (2D), or surface 
(3D). When the frequency of the two waves differs the shock moves with the average phase velocity, provided there 
are no phase slips in ID, or its equivalent in higher dimensions (creation of vortex pairs in 2D and inflation of vortex 
loops in 3D), for review see Bohr et al (1998). A stationary shock is formed most easily when a wave impinges on 
an absorbing boundary. Here we will consider the general situation where a plane wave is perturbed by a stationary, 
localized disturbance. This concept will be particularly useful in the context of interaction of defects. 

SufBciently upstream (i.e. against the group velocity) from the disturbance the perturbation will be small and we 
can linearize around the plane- wave solution. In general one obtains exponential behavior with exponents p = ik 
calculated from the dispersion relation (19) with A = 0. This gives 

p{(l + 62)(4Q2 +p-^)p _ 2F2[(1 + bc)p - 2(6 - c)Q]} = . (25) 

After separating out the translational mode p = one is left with a cubic polynomial. To discuss the roots pi, p2, P3 
choose the group velocity V^n = 2(6 — c)Q > (otherwise all signs must be reversed). For V^y = one has 
Pi = —pz < 0,p2 = in the Eckhaus stable range. For Vg\\ > one finds that p2 increases with increasing \Q\. 
Before \Q\ = 1 is reached p2 and ps collide and become complex conjugate. 

The existence of roots with positive real part is an indication of screening of disturbances in the upstream direction 
because one needs the solutions that grow exponentially to match to the disturbance. The screening length is given 
by {Rep2)~^, the root with the smaller real part. So the screening is in general exponential. The length diverges for 
^ and then one has a crossover to power law. 

In the focusing quadrants the root collision always occurs before the Eckhaus instability is reached. (It fails to hold 
only in part of the defocusing quadrant away from the origin and restricted to > 6^.) We will refer to the situation 
where p2,3 are real to the monotonic case and otherwise to the oscillatory case, because this characterizes the nature 
of the asymptotic interaction of sources that emit waves, see Sec. Ill and IV. In Fig. 2 the transition from monotonic 
to oscillatory behavior for the waves emitted by standing hole solutions is also shown (curve MOH) . A generalization 
to disturbances that move with velocity v is straightforward by replacing in Eq. (19) the growth rate A by —vp. 

G. Phase equations 

The global phase invariance of the CGLe leads to the fact that, starting from a coherent state (in particular a 
plane wave), one expects solutions where the free phase becomes a slowly varying function, and one can construct 
appropriate equations for these solutions. (In the case of localized structures one then also has to allow for a variation 
of the velocity.) In the mathematical literature these phase equations are sometimes called "modulated modulation 
equations" . There are several ways to proceed technically in their derivation, but the simplest is to first establish the 
linear part of the equation by a standard linear analysis and then construct nonlinearities by separate reasoning. It 
is helpful to include symmetry considerations to exclude terms from the beginning on. 
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By starting from a plane wave state (17) with wave number Q in the x direction one may perform a gradient 
expansion of </> leading in ID to 

dt(f> = Vg\\d^(l) + D2\\dl(t> - ^g\\dl<p - D^dt(t> + h.o.t. (26) 

Note the absence of a term proportional to which would violate translation invariance. The prefactors of the gradient 
terms have to be chosen according to (21) in order to reproduce the linear stability properties of the plane wave (with 
the wavenumber Qk replaced by Q). 

Nonlinearity can be included in Eq. (26) by substituting in the coefficients of Eq. (21) Q — > Q + dx<i>- Expanding 
Vg||((5) and D-2\\{Q) one generates the leading nonlinear terms, so that eq. (26) goes over into (Janiaud et al, 1992) 

{dt - Fg||9.)0 = D^\\dl<P - D^d^^cj) - 5i(5,</.)2 _ ^^^^q3^ _ g^Q^^Ql^ + h.o.t. (27) 

with the linear parameters from (21) and 



51 = -T^dqVgW = {h-c) 



1 

2' 

c, n 4go(l+c2) 
52 = -5qL»2|| = i-Ql (28) 

Because of translation invariance the lowest relevant nonlinearities are {d^(f>)'^ and dx(f>d^(f) [Kuaramoto, 1984]. 

In the stationary case (c = 6) one has V^|| = = gi = 0. The remaining nonlinearity dx4>dl(j) does not saturate 
the linear instability and one recovers the results for the nonlinear Eckhaus instability (Kramer and Zimmermann, 
1985). For b ^ c and Q = the dominant nonlinearity is {d^cf))^ which does saturate the linear instability. The 
resulting equation with flg^^ = g2 = 

dt(t>-D2iidl<j> + D^dt<p + gi{dx^f =0 (29) 

describes the (supercritical) bifurcation at the Benjamin-Feir-Newell instability and is known as the Kuramoto- 
Sivashinsky equation. It has stationary periodic solutions, which are stable in a small wavenumber range 0.77j3o < 
p < 0.837po7 where po = \fD^ilDl (Frisch et al. 1986, Nepomnyashchu, 1995b). It also has spatio-temporally chaotic 
sohitions, which are actually the relevant attractors, and it maybe represents the simplest partial differential equation 
displaying this phenomenon. Since it is a rather general equation for long-wavelength instabilities, independent of the 
present context of the CGLe, it has attracted much attention in recent years (see e.g. Bohr et al (1998)). 

In the general case Q ^ both nonlinearities are important (and also the term proportional to f2g). Then the 
bifurcation to modulated waves, which are represented by periodic solutions of (26), can be either forward or backward 
depending on the values of b, c. In fact the Eckhaus bifurcation always becomes supercritical before the Benjamin- 
Fcir-Newell instability is reached. Modulated solutions can be found analytically in the limit D2\\ — > 0, where Eq. 
(26) reduces to 

dt(t> - yg\\dx<t> + dl'P + {dx(t>f + 0(e) = 0. (30) 

with e = DiD2/\^g\ <C 1. This limit is complementary to that leading to the Kuramoto-Sivashinsky equation and 
one is in fact left with a perturbed Korteweg de Vries equation. This equation is (in ID) completely integrable and the 
perturbed case was studied by Janiaud et al (1992) and by Bar and Nepomnyashchii (1995). They showed that a fairly 
broad band of the periodic solutions (much larger than in the Kuramoto-Sivashinsky equation equation) persists the 
perturbation of the Korteweg de Vries equation and is stable. Outside this limit, and in particular in the Kuramoto- 
Sivashinsky equation regime, one has in extended systems spatio-temporally chaotic solutions. It is believed that this 
chaotic state is a representation of the phase chaos observed in the CGLe, although this view has been challenged (see 
Sec. III.D). Sakaguchi (1992) introduced higher-order nonlinearities into the Kuramoto-Sivashinsky equation, which 
allowed to capture the analog of the transition (or cross over) to amplitude chaos manifesting itself by finite-time 
singularities in dxip- 

Under some conditions, namely when D2\\ is positive and Q is small, the first nonlinearity in Eq. (27) is dominant 
and then it reduces to Burgers equation 

dt(t)-D2\\dl4> + gi{dx4>f = Q (31) 

which is completely integrable within the space of functions that do not cross zero, because it can be linearized by a 
Hopf-Cole transformation (j) = (Z)2||/fi(i) log W^. The equation allows to describe analytically sink solutions (shocks), 
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(sec e.g., Kuramoto and Tsuzuki (1976), Malomed (1983) and is (in 2D) useful for the qualitative understanding of 
the interaction of spiral waves in the limit 6, c ^ (see Biktashev (1989), Aranson et al (1991)). 

Clearly generalization of the phase equations to higher dimensions is possible. Using the full form of Eq. (21) one 
can easily generalize the phase equation to 2D (for details sec also Kuramoto (1984); Lega (1991)). The equation 
is most useful in the range where D^w goes through zero and becomes negative. Also, they can be derived outside 
the range of applicability of the CGLe (far away from threshold) and for more general nonlinear long-wavelength 
phenomena, where (j) does not represent the phase of a periodic function (see e.g. Bar and Nepomnyashchii (1995)). 

H. Topological defects 

Zeros of the complex field A result in singularity of the phase 6 = arg^. In 2D point singularities correspond 
to quantized vortices with so-called topological charge n = l/(27r) V^dl, where L is a contour encircling the zero 
of A. Although for & 7^ c they represent wave emitting spirals, they are analogous to vortices in superconductors 
and superfluids and represent topological defects because a small variation of the field will not eliminate the phase 
circulation condition. Clearly, vortices with topological charge n = ±1 are topologically stable. The vortices with 
multiple topological charge can be split into single-charged vortices ^. 2D point defects become line defects in 3D. 
Then, one can close such a line to a loop, which can shrink to zero. Some definitions of topological defects from the 
point of view of energy versus topology are given by Pismen (1999). 

Topological arguments do not guarantee the existence of a stable, coherent solution of the field equations. In 
particular, the stability of the topological defect depends on the background state it is embedded in. For example, spiral 
waves are stable in a certain parameter range, where they select the background state (see Sec. IV). Simultaneously 
charged sinks coexist with spirals and play a passive role. Also, defects can become unstable against spontaneous 
acceleration of their cores. 

I. Effects of boundaries 

Boundaries may play an important role in nonequilibrium systems. Even in large systems the boundary may provide 
restriction or even selection of the wavenumber. We will not discuss this topic (see, e.g. Cross and Hohenberg, 1993) 
and consider situations where such effects are not important. 

III. DYNAMICS IN ID 

In this Section we will consider the properties of various ID solutions of the CGLe. We will discuss the stability 
and interaction of coherent structures and the transition to and characterization of spatio-temporal chaos. 

A. Classification of coherent structures, counting arguments 

Coherent structures introduced in the previous chapter can be characterized using simple counting arguments put 
forward by van Saarloos and Hohenberg (1992). It should be mentioned, however, that the counting arguments 
cannot account for all the circumstances, for example, hidden symmetries, and may fail for certain class of solutions, 
see below. ID coherent structures can be written in the form 

A{x, t) = e-*'^'=*+*'?^(^-''*)a(.T - vt), (32) 

with the real functions </>, a satisfying a set of 3 ordinary differential equations 

ax = s; 

Sx = fA^o — 7^"^[(1 + bujk)a + v{s + b^a) — (1 + bc)a^]; 

tpx = —2sip/a + 7~^[6 — LUk+ v{sb/a — ip) — {b — c)a^]. (33) 



^In small samples vortices with multiple charge can be dynamically stable, see for detail Geim et al (1998); Deo et al (1997). 
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Abbreviations used are •j = 1 + b"^ , = da/dx,Sx = ds/dx and tf) = d(j)/dx. These ordinary differential equations 
constitute a dynamical system with three degrees of freedom. 

The similarity transform of Sec. II. C can be adapted to these equations by absorbing the external wavenumber in 
the phase (j), which leads to Wfe = a; — Qv, lo'^,=w' — Q'v', and 

, . (34) 

The other relations (14-16) remain. 

The counting arguments allow to establish necessary conditions for the existence of localized coherent structures, 
which correspond to homo- or heteroclinic orbits of Eqs. (33). Consider, e.g. the trajectory of Eqs. (33) flowing 
from fixed point to fixed point L. If has n^y unstable directions, there are n^r — 1 free parameters characterizing 
the flow on the njv-dimensional subspace spanned by the unstable eigenvectors. Together with the parameters 
and V this will yield njv + 1 free parameters. If L has unstable directions, the requirement that the trajectory 
should come in orthogonal to these yields conditions. The multiplicity of this type of trajectory will, therefore, be 
n = Un — TiL + I, and, depending on n, it will give either a n-parameter family (n > 1), a discrete set of structures 
(n = 0) or no structure (n < 0). In addition one may have symmetry arguments that reduce the number of conditions. 

The asymptotic states can either correspond to non-zero steady states (plane waves), or to the trivial state a — 0. 
Accordingly, the localized coherent structures can be classifled as pulses, fronts, domain boundaries, and homoclons. A 
pulse corresponds to the homoclinic orbit connecting to the trivial state a = 0. They come in discrete sets (Hohenberg 
and van Saarloos, 1992). For a full description of the pulses sec Akhmcdicv et al (1995,1996,2001), Afanasjcv et al 
(1996). Fronts are heteroclinic orbits connecting on one side to a plane- wave state and to the unstable trivial state 
on the other side. They come in a continuous family, but sufficiently rapidly decaying initial conditions evolve into a 
"selected" front that moves at velocity v* = 2\/\ + 6^ and generates a plane wave with wave number Q* ~ b/^/1 + 6^ 
("linear selection"). For a discussion, see Hohenberg and van Saarloos (1992), Cross and Hohenberg (1993). 

Domain boundaries are heteroclinic orbits connecting two different plane-wave states. The domain boundaries can 
be active (sources, often called holes) or passive (sinks or shocks) depending on whether in the co-moving frame the 
group velocity is directed outward or inward, respectively. They will be discussed in the next subsection. Homoclons 
(or homoclinic holes, phasons) connect to the same plane-wave state on both sides. They are embedded in solutions 
representing periodic arrangements of homoclons, which represent quasiperiodic solutions satisfying ansatz (32) and 
correspond to closed orbits of Eqs. (33). They will be discussed in Subsec. III.C.2. There also exist chaotic solutions 
of Eqs. (33), which correspond to nonperiodic arrangement of holes and shocks, or homoclons, see Subsecs.III.B.2, 
III.C.2. 

B. Sinks and sources, Nozaki-Bekki hole solutions 

Sinks (shocks) conform with the counting arguments of Hohenberg and van Saarloos (1992): there is a two-parameter 
family of sinks. However, there are no exact analytic expressions for the sink solution connecting two traveling waves 
with arbitrary wavenumbers. The exact sink solution found by Nozaki and Bekki (1984) corresponds to some special 
choice of the wavenumbers and is not therefore typical. For small difference in wavenumber a phase description of 
sinks is possible (Kuramoto, 1984). 

However, the counting arguments (which arc, strictly speaking, neither necessary nor sufficient, and cannot account 
for speciflc circumstances, such as a hidden symmetry) fail for the source solution. According to the counting 
arguments there should be only a discrete number, including in particular the symmetric standing hole solution that 
has a zero at the center and emits plane waves of a definite wavenumber (for given 6, c). However, the standing hole is 
embedded in a continuous family of analytic moving sources, the Nozaki-Bekki hole solutions. They are characterized 
by a localized dip in \A\ that moves with constant speed v and emits plane waves with wave numbers q\ qi- 

This is a special feature of the cubic CGLe as demonstrated by the discovery that the moving holes do not survive a 
generic perturbation of the CGLe, e.g. a small quintic term (see below). Thus they are not structurally stable (Popp 
et al, 1993, 1995, Stiller et al, 1995a,b). With the perturbation the holes either accelerate or decelerate depending on 
the sign of the perturbation and other solutions appear, see also Doleman (1995). Clearly, the ID CGLe possesses a 
"hidden symmetry" and has retained some remnant of integrability from the NLSe. Apparently this non-genericity 
has no consequences for other coherent structures. 

The Nozaki-Bekki hole solutions are of the form (Nozaki and Bekki, 1984) 



y, = b' + p'i- 



, 1 b - b' 

uik-b+ V 

2 7 
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A^^ = [Bd(;(pv{KO + Av] X exp[i(^„(KC) + iav - iflt] (35) 
where ¥'v(kC) = liicosh(K;^) 

and ~ X — vt, V is the velocity of the hole, and (ji/2 are the asymptotic wavenumbers. Symbols with a "hat" denote 
real constants depending only on b and c, e.g. d = 1/(2^— c)) . The frequency Q and are linear functions of v^. 
The emitted plane waves have wavenumbers 

gi/2 = ±(3 + a (36) 

where (3 = k/k and a = va. One easily derives the relation 

V={uj{q2)-uj{qi))/{q2-qi) (37) 

where w(g) is the dispersion relation for the plane waves. This relation can be interpreted as phase conservation and is 
valid also for more general equations possessing phase invariance. For the cubic CGLe it reduces to tj = (6 — c) ((71 + (72), 
i.e. the hole moves simply with the mean of the group velocities of the asymptotic plane waves. The exact relations 
between the parameters can be derived by inserting ansatz (35) into the CGLe. The resulting algebraic equations 
(8 equations for 8 parameters) turn out to be not independent, yielding the one-parameter family, becomes zero 
at a maximal velocity -iiVmax and here the hole solution merges with a plane wave with wavenumber qi — q2 = 
Vmax / {'2'(b — c)) ■ In a large part of the b,c plane these bifurcations occur in the range of stable plane waves. This is 
the case outside a strip around the line b = c given by |6 — c| > i5 with S varying from 1/2 for large values of |&|, |c| 
to 0.55 for small values. The region extends almost to the Benjamin- Feir-Newell curves. The bifurcation cannot be 
captured by linear perturbation around a plane wave, since the asymptotic wavenumbers of holes differ. Presumably 
the bifurcation can be captured by the phase equation (27). For b c only the standing hole survives (the velocities 
of the moving holes diverge in this limit). 

The interest in the hole solutions comes particularly from the fact that they are dynamically stable in some range. 
The stability has first been investigated by Sakaguchi (1991), in direct simulations of the CGLe. The stability problem 
was then studied by Chate and Manneville (1992), numerically for u = 0,u = 0.2 and by Sasa and Iwamoto (1992), 
semi-analytically for v = and by Popp et al. (1995), essentially analytically. As a result, hole solutions were found 
to be stable in a narrow region of the b-c plane which is shown in Fig. 2 for the standing hole {v = 0) (upper shaded 
region). From below, the region is bounded by the border of (absolute) stability of the emitted plane waves with 
wavenumber q{b,c) (see Eq.(35)) corresponding to the continuous spectrum of the linearized problem, see curve AH 
in Fig. 2. From the other side, the stable range is bounded by the instability of the core with respect to localized 
eigenmodes corresponding to a discrete spectrum of see curve CI in Fig.2. For c — > 00, one has 

bci for 6 > (38) 

bci^ — ^ + 0(l/c) for6<0 
v2 

The result could be reproduced fully analytically by perturbing around the NLSe limit, where the Nozaki-Bekki holes 
emerge as a subclass of the 3-parameter family of dark solitons and is supported by detailed numerical simulations 
and shows excellent agreement (Stiller et al, 1995a). It differs from that obtained by Lega and Fauve (1997) and 
by Kapitula and Rubin (2000). Lega and Fauve (1997), Lega (2000) claim a larger stability domain which is in 
disagreement with simulations (Stiller et al, 1995a). 

The core instability turns out to be connected with a stationary bifurcation where the destabilizing mode passes 
through the neutral mode related to translations of the hole. This degeneracy is specific to the cubic CGLe and is 
thus structurally unstable (see below). When going through the stability limit the standing hole transforms into a 
moving one. Indeed, the cores of moving holes were found to be more stable than those of the standing ones (Chate 
and Manneville, 1992). 

1. Destruction of Nozaki-Bekki lioles by small perturbations 
Consider the following "perturbed cubic CGLe" 

dtA = [1 + (1 + ib)dl - (1 + ic)\Af + d\A\^] A (39) 



14 



where a quintic term with a smaU complex prefactor d = d' + id" {\d\ ^ 1) is included, which will be treated as a 
perturbation. There are, of course, other corrections to the cubic CGLe, but their "perturbative effect" is expected 
to be similar. 

Simulations with small but finite d, show that stable, moving holes are in general either accelerated and eventually 
destroyed or slowed down and stopped to the standing hole solution depending on the phase of d (Popp et al 1993, 
1995; Stiller et al 1995a,b). In particular, for real d = d' one has 



dtv 



V 



^O^^d'^O. (40) 



One finds that the relations (37) connecting the core velocity and the emitted wavenumbers are (almost) satisfied 
at each instant during the acceleration process. The acceleration thus occurs approximately along the Nozaki-Bekki 
hole family and it can be described by taking v = v{t) as a slowly varying variable while other degrees of freedom 
follow adiabatically. The semianalytic matching-perturbation approach gives the reduced acceleration, which agrees 
with the simulations, see Fig. 3. 

Of special interest is the case where the core-stability line is crossed while d (which corresponds to the typical 
situation with higher-order corrections to the CGLe). Then the two modes which cause the acceleration instability 
and the core instability (which is stationary for d = 0) are coiiplcd which in the decelerating case leads to a Hopf 
bifurcation. As a result, slightly above the (supercritical) bifurcation, one has solutions with oscillating hole cores 
(Popp et al 1993), see Fig. 4d. The normal form for this bifurcation - valid for small \d\,v,u - is 

u = {X - sv'^)u + div (41) 
V = IJ,u + d2V , (42) 

which is easy to analyze. Here u and v are the amplitudes of the core-instability and acceleration-instability modes 
respectively, s and /x are of order dP while di^2 must be of order d^, since in the absence of a perturbation holes with 
nonzero velocity exist. A can be identified with the growth rate of the core-instability mode at d = 0. The nonlinear 
term in Eq.(41) takes care of the fact that moving holes are more stable than standing ones and at the same time 
saturates the instability (s > 0). 

Far away from the core-instability threshold where A is strongly negative u can be eliminated adiabatically from 
Eqs.(41) which ior v ^ yields 

^ = (^2-^1^)^; (43) 

The term in brackets can be identified with the growth rate of the acceleration instability. The parameters of Eqs. 
(43) were calculated fully analytically by Stiller et al (1995a) for 6, c — > oo, 6c > 0. 

2. Arrangements of holes and shocks 

When there are more than one hole they have to be separated by shocks. The problem of several holes is analogous 
to that of interacting conservative particles, which is difficult to handle numerically. The solutions resulting from a 
periodic arrangement of holes and shocks are actually special cases of the quasiperiodic solutions - "homoclons" in 
the limit of large periodicity - to be discussed in the next subsection. In the situation discussed here the solutions 
depend sensitively on perturbations of the CGLe 

Such states are frequently observed in simulations with periodic boundary conditions (see e.g. Manneville and 
Chate (1992), Popp et al (1993), Stiller et al (1995)a,b). Figure 5 shows the modulus \A\ = \A{x)\ of a typical solution 
found in a simulation. As shown in Figure 4 one finds uniform as well as (almost) harmonic and strongly anharmonic 
oscillating hole velocities. (Slightly) beyond the core instability line the direction of the velocity is changed in the 
oscillations (see Figure 4 d). The solutions are seen to be very sensitive to d-perturbations of the cubic CGLe. 

The uniformly moving solutions can be well understood from the results of the last subsections. First, they are not 
expected to exist (and we indeed could not observe them in simulations) in the range of monotonic interaction, below 



*For d = one expects to observe conservative dynamics for holes and shocks, similar to interacting particles for zero friction. 
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curve MOH in Fig. 2, ^ since here the asymptotic hole-shock interaction is always attractive. In the oscillatory range, 
and away from the core instability line uniformly moving periodically modulated solutions can then be identified as 
fixed points of a first-order ordinary differential equation for the hole velocity v which can be derived by the matching- 
perturbation method (Popp et al, 1995). In addition, solutions with oscillating hole velocities were found coexisting 
(stably) with uniformly moving solutions (see Figures 4b, 4c). They can (in a first approximation) be identified as 
stable limit cycles of a two-dimensional dynamical system with the hole velocity v and the hole-shock distance L as 
active variables. These oscillations are of a different nature than those found slightly beyond the core instability line 
(Fig. 4d), see above. 

By an analysis of the perturbed equations (33) Doelman (1995) has shown that the quintic perturbations create 
large families of traveling localized structures which do not exist in the cubic case. 

3. Connection with experiments 

Transient hole-type solutions were observed experimentally by Lega et al (1992), and Flesselles et al (1994) in the 
(secondary) oscillatory instability in Rayleigh-Benard convection in an annular geometry. Here one is in a parameter 
range where holes are unstable in the cubic CGLe, so that small perturbations are irrelevant. Long time stable sta- 
tionary holes ("Id spirals") were observed in a quasi-ld chemical reaction system (CIMA reaction) undergoing a Hopf 
bifurcation by Perraud et al (1993). The experiments were performed in the vicinity of the cross-over (codimension-2 
point) from the (spatially homogeneous) Hopf bifurcation to the spatially periodic, stationary Turing instability. Sim- 
ulations of a reaction-diffusion system (Brusselator) with appropriately chosen paranictc;rs exhibited the hole solutions 
(and in addition more complicated localized solutions with the Turing pattern appearing in the core region). Strong 
experimental evidence for Nozaki-Bekki holes in hydrothermal nonlinear waves is given by Burguette et al (1999). 

Finally we mention experiments by Leweke and Provensal (1994), where the CGLe is used to describe results of 
open-flow experiments on the transitions in the wake of a bluff body in an annular geometry. Here the sensitive 
parameter range is reached and in the observed amplitude turbulent states holes should play an important role. 

C. Other coherent structures 

Coherent states (32) with periodic hmctions a and ^ = d(j)/dx (same period) have been of particular interest. In 
general, when the spatial average of tp is nonzero, the complex amplitude A is quasiperiodic, i.e. A can be written 
in the form (12) with B a periodic function of ^ = .x — vt. The associated wavcnumbcr p will be called the "inner 
wavenumber" , in contrast to the outer wavcnumbcr Q, which is equal to the spatial average of ip. These quasiperiodic 
solutions bifurcate from the traveling waves (17) in the Eckhaus unstable range Q|, < < 1 at the neutrally stable 
positions obtained from Eq. (19) with ReX = and k replaced by ]3(|p| < 2). In fact, the long-wave Eckhaus instability 
is signalizes in particular by bifurcations of the solitary (or homoclinic) limit solution p ^ 0. In this limit the velocity 
V is at the bifurcation equal to the group velocity Vg, see Eq. (21). 

Well away from the Benjamin-Feir-Newell instability (on the stable side) the bifurcation is subcritical (as for the 
GLe). However, it becomes supercritical before the Benjamin-Feir-Newcll line is reached (Janiaud et al (1992)) and 
remains so in the unstable range. The bifurcation is captured analytically most easily in the long-wave limit by phase 
equations, see Subsec. II. G. The supercritical nature of the bifurcation allows to understand the existence of "phase 
chaos" that is found when the Benjamin-Feir-Newell line is crossed (see below). 

1. The GLe and NLSe 

For the GLe a full local and global bifurcation analysis is possible and the quasi-periodic solutions can be expressed 
in terms of elliptic functions (Kramer and Zimmermann, 1985; Tuckerman and Barkley, 1990). Indeed, for b = c = 
Eqs. (33) lead to the second-order system 

axx = daU where U = ^[-a^ + ^a^ - h^a^ (44) 



^The boundary between monotonic and oscillatory interaction depends on the hole velocity (see Subsec. II. F). 
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with h = a^tp an integration constant. This allows to invoke the mechanical analog of a point particle (position a, 
time X ) moving in the potential —U{a). On sees that for \h\ < He = •\/4/27 « 0.385 the potential — J7 has extrema 

at ai,2 with 012^! — 2 = ^ corresponding to the plane-wave solutions (17) with Fi^2 = ^1,2, \Qi,2\ = ^1 — 2- 

The maximum corresponds to \Qi\ < Qe '■= l/-\/3. It is stable since it correspond to a minimum of V. The solution 
IQ2I > is unstable. In this way the Eckhans instability is recovered. The bifurcations of the quasi-periodic 
solutions are subcritical and the solutions are all unstable. They represent the saddle points separating (stable) 
periodic solutions of different wavenumber, and thus characterize the barriers against wavelength-changing processes 
involving phase slips. 

As was already mentioned in the Sec. II, all stationary solutions and their stability properties of the ID (defocusing) 
NLSe coincides with those of the GLe. In addition, the NLSe has more classes of coherent structures due to the 
additional Galilean and scaling invariance absent in GLe. We will not touch this question since the ID NLSe is a 
fully integrable system (Zakharov and Shabat, 1971) that has been studied in great detail (see, e.g. Proceedings of 
Conference on The Nonlinear Schrodinger Equation (1994)). 



2. The CGLe 

For the CGLe the local bifurcation analysis has been first carried out in the limit of small p using the phase equation 
(27) (Janiaud et al 1992) and subsequently for arbitrary p (Hager and Kramer, 1996). It shows that the Eckhaus 
instability becomes supercritical slightly before the Benjamin-Feir-Newell curve is reached. The bifurcation has in the 
long-wave limit a rather intricate structure, since the limits p ^ and Q Qe do not interchange (this can already 
be seen in the GLe). When taking first the solitonic limit p ^ (while Q > Qe) the Eckhaus instability becomes 
supercritical slightly later than in the other ("standard") case, which corresponds to harmonic bifurcating solutions. 
These features are captured nicely by the phase equation description. 

Recently Brusch et al (2000) carried out a systematic numerical bifurcation analysis based on Eqs. (33) for the 
case Q = 0, where A is periodic (nevertheless we will usually refer to these solutions as quasi-periodic solutions). 
It shows that the supercritically bifurcating quasi-periodic branch {shallow quasi-periodic solutions) terminates in a 
saddle-node bifurcation and merges there with an "upper" branch {deep quasi-periodic solutions), see Fig. 6. For 
p 7^ the bifurcating solution has velocity v = (in Eqs. (33) the bifurcation is of Hopf type). It is followed by a 
drift-pitchfork bifurcation generating a nonzero velocity. The separation between the two bifurcations tends to zero for 
vanishing p, so that the solitary solutions develop a drift from the beginning on (at second order in the amplitude). ^ 
These features can be reproduced by the Kuramoto-Sivashinsky equation, which only exhibits the shallow homoclons. 
By adding higher-order (nonlinear) terms (Sakaguchi, 1990) the other branch can be generated. Brusch et al (2000) 
give evidence that the existence of the two branches provides a mechanism for the stabilization of phase chaos (see 
below). 

Deep quasi-periodic solutions (in the solitary limit) were first studied by van Hecke (1998) and then by van Hecke 
and Howard (2001) in the context of spatio-temporal chaos in the intermittent regime (see next subsection). Following 
this author we have adopted the name "homoclons" for the localized objects. 

The stability properties of quasi-periodic solutions were also analyzed by Brusch et al, 2000. Both branches of 
quasi-periodic solutions have neutral modes corresponding to translation and phase symmetries. The eigenvalue 
associated with the saddle-node bifurcation is positive for deep homoclons and negative for shallow ones. Apart from 
these three purely real eigenvalues, the spectrum consists of complex conjugated pairs. 

For not too small wavenumber p all the eigenvalues of the shallow branch are stable within a system of length 
L = 2'jt/p, but when L increases, the quasi-periodic solutions may become unstable with respect finite-wavelength 
instabilities (near the bifurcation they certainly do). 

Shallow quasi-periodic solutions can be observed in simulations of the CGLe (away from the bifurcation from the 
plane waves), in particular for Q ^ (Janiaud et al, 1992; Montague, et al 1996, 1997; Torcini, 1996; Torcini, 
Praunkron, and Grassberger, 1997). 

Also more complex coherent structures, corresponding to nonperiodic arrangements of shallow homoclons were 
found numerically. They can be understood from the fact that the interaction between shallow homoclons can be of 



for Q ^ there is drift already after the first bifurcation, so the drift-pitchfork gets unfolded (Brusch, private communication). 
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oscillatory nature depending on parameters. This is suggested by the fact that in the supercritical regime the spatial 
exponents introduced in Sec. II. F are complex. Clearly chaotic solutions may be expected to exist in Eqs. (33). 

Recently, quasi-periodic solutions have been observed experimentally in the form of modulations of spiral waves 
excited by core meandering (see Sec. IV.B.2,IV.H.2) in an oscillatory reaction-diffusion system by Zhou and Ouyang 
(2000). 

D. Spatio-temporal chaos 

1. Phase chaos and the transition to defect chaos 

When crossing the Benjamin-Feir-NewcU line with initial conditions A 1 (actually any nonzero, spatially constant 
A is equivalent) one first encounters phase chaos, which persists approximately in the lower dashed region between 
Benjamin- Feir-Newell and absolute instability curves in Fig. 2. In this spatio-temporally chaotic state \A\ remains 
saturated (typically above about 0.7 — 0.9) so there is long-range phase coherence and the global phase difference 
(here equal to 0) is conserved. In this parameter range one also has stable, periodic solutions, but they obviously have 
a small domain of attraction which is not reached by typical initial conditions. Beyond this range phase slips occur 
and a state with a nonzero (average) rate of phase slips is established (defect or amplitude chaos). Since in phase 
chaos only the phase is dynamically active it can be described by phase equations (in the case of zero global phase 
the the Kuramoto-Sivashinsky equation, otherwise Eq. (27)). 

In the CGLe, phase chaos and the transition to defect chaos (see below), was first studied by Sakaguchi (1990) and 
subsequently systematically in a large parameter range by Shraiman et al (1992) and for selected parameters and in 
particular very large systems by Egolf and Greenside (1995). One of the interests driving these studies was the question 
if in phase chaos the rate of phase slips is really zero, so it could represent a separate phase (in the thermodynamic 
sense), or if the rate is only very small. The recent studies by Brusch et al (2000) on the quasi-periodic solutions with 
small modulation wavenumber p (see subsec. III.C.2) point to a mechanism that prevents phase slips in a well-defined 
parameter range. 

In fact, phase chaos "lives" on a function space spanned by the forwardly bifurcating branches of the quasi- 
periodic solutions. Thus, snapshots of phase chaos can be characterized (roughly) as a disordered array of shallow 
homoclons (or "disordered quasi-periodic solutions"), see Fig. 7. The dynamics can be described in terms of birth and 
death processes of homoclonal units. When these branches are terminated by the saddle-node bifurcation, phase-slip 
processes are bound to occur. Brusch et al (2000) have substantiated this concept by extensive numerical tests with 
system sizes ranging from L = 100 to L = 5000 and integration times up to 5 x 10^. For a given L (not too small) 
phase slips occurred only past the saddle-node for the quasi-periodic solutions with inner wavenumber p = 2tt/L. 
Tracking of (rare) phase-slip events corroborated the picture. Thus, the authors conjectured that the saddle-node line 
for p ^ Q provides a strict lower boundary for the transition from phase to defect chaos. 

One expects existence of a continuous family of two types of phase chaos with different background wavenumbers 
Q. For Q ^ the state should arise when crossing the Eckhaus boundary for the plane wave with wavenumber Q, 
before the Benjamin-Feir-Newell limit. Such states were studied by Montague et al (1996,1997) and Torcini (1996) 
and Torcini et al (1997). It was found that the parameter range exhibiting phase chaos decreases with increasing 
Q. Thus, at fixed parameters b, c, one has a band of phase chaotic states which is bounded from above by some 
Qmax(b, c). Approaching the limit of phase chaos, Qmax decreases smoothly to zero, i.e. the phase chaotic state with 
Q = is the last to loose stability, similar to the situation for the plane waves. As mentioned before, for not too small 
Q, one finds coexistence with quasi-periodic solutions and spatially disordered coherent states. 

2. Defect chaos 

The transition (or crossover) between phase and defect chaos is reversible only for |6| larger than about 1.9 (near 
the lower edge of Fig. 2) [Shraiman et al, (1992)]. There, when approaching the transition from the side of defect 
chaos, the phase slip rate goes smoothly to zero. On the other hand, for \b\ < 1.9, there is a region where phase and 
defect chaos coexist ("bichaos", lower dashed region to the right of line DC), and in fact defect chaos even persists 
into the Benjamin-Feir-Newell stable range (the limit towards small c is approximated by the dashed line DC in Fig. 
2) (Shraiman et al, 1992; Chate, 1993; Chate and Mannneville, 1994). There the defect chaos takes on the form of 
spatio-temporal intermittency. 
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A qualitative understanding of the parameter region where one has bistable defect chaos can be obtained as follows: 
In this region, starting from the ("saturated") state with \A\ « 1 and forcing a large excursion of |A| leading to a 
phase slip, makes it easier for other phase slips to follow. This memory effect is suppressed with increasing |6| as can 
be seen from the increase of the velocity v* = 2\/l + b'^ with \b\ of fronts that tend to restore a saturated state (see 
Sec. III. A). 

Prom above (towards small values of \b\) the chaotic state joins up with the region of stable Nozaki-Bekki holes. 
The characteristics of these holes is influenced by small perturbations of the CGLe, and this in turn affects the precise 
boundary of of spatio-temporal chaos, see Subsec. III.D.4. 

3. The intermittency regime 

This regime where defect chaos coexists with stable plane waves has been studied numerically in detail by Chate 
(1993, 1994) and he has pointed out the relation with spatio-temporal intermittency. There, typical states consist of 
patc;hes of plane waves, separated by various localized structures characterized by a depression of \A\. The localized 
structures can apparently be divided into two groups depending on the wavenumbers qi and Qr of the asymptotic waves 
they connect. On one hand one has slowly moving structures which can be related to Nozaki-Bekki holes, which in 
this regime are either core stable or have long lifetime. They become typical as the range of stable Nozaki-Bekki holes 
is approached (Chate, 1993, 1994). 

However, in a larger range the dominant local structures have velocities and asymptotic wavenumbers that are 
incompatible with the Nozaki-Bekki holes and which, according to van Hecke (1998), must be associated with the 
deep homoclons. Space-time plots of the amplitude |^|, the phase (arg(^)), and the local wavenumber q:=dx arg(A) 
in such a regime are shown in Fig. 8 (the plot is a zoom in of a larger simulation as shown in Fig. 9b). The 
wavenumbers of the laminar patches arc quite close to zero, while the cores of the local structures are characterized 
by a sharp phase-gradient (peak in q) and dip of |^|). The holes propagate with a speed of 0.95 ± 0.1 and either their 
phase gradient spreads out and the hole decays, or the phase-gradient steepens and the hole evolves to a phase slip. 
The phase slip then sends out one (or maybe two) new localized objects which repeat the process. Van Hecke (1998) 
provides evidence that the dynamics evolves for much of the time around the one-dimensional unstable manifold of 
homoclons with background wavenumber near zero. They provide the saddle-points that separate dynamical processes 
which lead to phase slips from those that lead back to the (laminar) plane-wave state. When the parameters h and c 
are quenched in the direction of the transition to plain waves, these zigzag motions of the holes becomes very rapid 
(Fig. 9a). 

Additional consideration of homoclons and their role in spatio-temporal chaos is presented by van Hecke and 
Howard, 2001. The simulations of the CGLe show that when an unstable hole invades a plane wave state, defects 
are nucleated in a regular, periodic fashion, and new holes can then be born from these defects. Relations between 
the holes and defects obtained from a detailed numerical study of these periodic states are incorporated into a simple 
analytic description of isolated "edge" holes, which are seen in Fig. 9. 

Recently Ipsen and van Hecke (2001) found in long-time simulations that in restricted parameter range composite 
zigzag patterns formed by periodically oscillating homoclinic holes represent the attractor. 

4. The boundary of defect chaos towards Nozaki-Bekki holes 

As the boundary of stability of Nozaki-Bekki holes is approached (curve AH in Fig. 2) one can increasingly observe 
Nozaki-Bekki hole-like structures that emit waves and thereby organize ( "laminarize" ) their neighborhood. In this 
regime the state becomes sensitive to small perturbations of the cubic CGLe, as discussed in Sec. III.B. 

It is found that stable hole solutions suppress spatio-temporal chaos and as a conseqiicnce, for a decelerating {d' < 
if d" = 0) perturbation, the (upper) boundary of spatio-temporal chaos is simply given by the stability boundary AS 
in Fig. 2 of the Nozaki-Bekki hole solutions, whereas for an accelerating perturbation spatio-temporal chaos is also 
observed further up (Popp et al 1993; Stiller et al 1995b). For d' < random initial conditions lead to an irregular 
grid of the standing holes, separated by shocks from each other, sec Fig. 10b. Such grids are the ID analog of the 
"vortex glass" state found in 2D (see Subsec. IV. G), see also Chate, 1993. 

Clearly, for accelerating d' > {] such a grid is unstable. In this case Nozaki-Bekki holes are created from random 
initial conditions, too. In their neighborhood they suppress small-scale variations of \A\ typical for amplitude chaos 
in the major part of the chaotic regime. However, since they are accelerated they only have finite life time, see 
Figure 10a. For parameters 6, c below about curve EH in Fig. 2 (the coincidence with the Eckhaus limit presumably 
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fortuitous) the destruction of the holes leads to creation of new holes and shocks, yielding a chaotic scenario of 
subsequent acceleration, destruction and creation processes. 

IV. DYNAMICS IN 2D 

A. Introduction 

The 2D CGLe has a variety of coherent structures. In addition to the quasi-lD solutions derived from the coherent 
structures in ID discussed above the CGLe possesses localized sources in 2D known as spiral waves. An isolated spiral 
solution is of the form 

A{r, e, t) = F{r) exp[i(-wi + raO + iP{r))] (45) 

where r, 9 are polar coordinates. The (nonzero) integer m is the topological charge, and u is the (rigid) rotation 
frequency of the spiral, F{r) > is the amplitude and xp{r) is the phase of the spiral, Q = drip for r ^ oo is the 
asymptotic wavenumber selected by the spiral (Hagan 1982), and ui — c + {b — c)Q'^ is the spiral frequency ^. Since 
spirals emit asymptotically plane waves (group velocity Vg = 2(6 — c)Q > 0), they are source solution. In addition 
to spirals there exist sinks which absorb waves. The form of sinks is determined by the configuration of surrounding 
sources. Sinks with topological charge appear as edge vortices (see below). Spiral solutions with m 7^ ±1 are unstable 
(Hagan, 1982). Single-charged spiral solutions are dynamically stable in certain regions of parameter space. The 
grey-coded image of a spiral and an edge vortex is shown in Fig. 11. For 6 — c — > 0, and also for 6, c — > 00, the 
asymptotic wavenumber of the spiral vanishes and the solution goes over into the well-known vortex solution of the 
GLe and NLSe, known in the context of superfluidity (Pitaevskii,1961, Gross 1963, Donnelly, 1991) and somewhat 
similar to that found in superconductivity in the London limit (Abrikosov, 1988, Blatter et al 1994). In periodic 
patterns spirals and vortices become dislocations. 

The asymptotic interaction is very different for the case 6 = c, where it it is long range decaying like with some 
corrections, whereas for 6 — c 7^ it is short range decaying exponentially (see below). Interaction manifests itself in 
a motion of each spiral. The resulting velocity can have a radial (along the line connecting the spiral cores) and a 
tangential component. 

Spirals may form regular lattices and/or disordered quasi-stationary structures called vortex glass or frozen state. 
When individual spirals become unstable (spiral break-up), the typical spatio-temporal behavior is chaotic. 

B. Spiral stability 

1. Outer stability 

For spirals to be stable the wavenumber of the asymptotic plane wave has to be in the absolutely stable range (see 
Sec. II. E). In order to find this stability limit we have to evaluate the condition (24) for the wavenumber Q emitted 
by the spiral, see Fig. 12. Existence of an absolutely stable spiral solution guarantees that small perturbations within 
the spiral will decay, but does not assure that it will evolve from generic initial conditions. For further discussion see 
Sec.IV.H.2. 

2. Core Instability 

The spiral core may become unstable in a parameter range where diffusive effects are weak compared to dispersion 
(large b limit). Then it is convenient to rewrite CGLe in the form 

dtA = A+{e + i)/:^A-{l + ic)\A\^A. (46) 



'^Spirals with zero topological charge m, — (targets) are unstable in the CGLe. However, inhomogeneities can stabihze the 
target. Hagan (1981) have found stable targets in the inhomogencous CGLe in the limit of small c. Stationary and breathing 
targets in the wide parameter range of the CGLe were studied by Hendrey et al (2000) 
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where e = l/\b\ and length has been rescaled by ^/e. This situation is typical in nonlinear optics (transversely extended 
lasers or passive nonlinear media). In this case a systematic derivation of the CGLe from the Maxwell-Bloch equations 
in the "good cavity limit" for positive detuning between the cavity resonance and the atomic line leads to very small 
values of £ > (CouUet et al 1989a; Oppo et al 1991; Newell and Moloney, 1992; Newell, 1994). Representative values 
are e ~ IQ-^ - IQ-^ (Coullet et. al, 1989a). 

For e = one has the Galilean invariance mentioned in Sec. II. C, (see Saarloos and Hohenberg, 1992, for the ID 
case) and then in addition to the stationary spiral there exists a family of spirals moving with arbitrary constant 
velocity v 

A{r, t) = F{r') exp i [-w't + 6 + ip{r') - (47) 

where r' = r + vt, w' = w — and the functions F, tp are those of Eq. (45) (this invariance holds for any stationary 
solution). For £ / the diffusion term ^ eAA destroys the family and in fact leads to ac- or deceleration of the 
spiral proportional to ev. The natural assumption is that one has deceleration so that the stationary spiral is stable 
(Coullet et. al, 1989a). In fact this is not the case. Stable spirals exist only above some critical value Ec- Below Sc 
the stationary spirals are unstable with respect to spontaneous acceleration (Aranson et al, 1994). 

For small values of e one may expect the solution (47) to be slightly perturbed and have a slowly varying velocity 
V which obeys an equation of motion of the form 9(V + eK\r = 0. Because of isotropy the elements of the tensor K 
must satisfy K^x = Kyy and K^y = —Ky^, so the relation can also be written as 

dtv + eKV = (48) 

where v = Vy + ivx and k = K^x — iK^y Since in general the friction constant k is complex the spiral core moves on 
a (logarithmic) spiral trajectory. 

The acceleration instability of the spiral core has a well-known counterpart in excitable media, where the spiral " tip" 
can perform a quasiperiodic motion leading to meandering (see, e.g. Barkley, 1994). The main difference between the 
two cases can be understood by considering the nonlinear extension of Eq. (48) dtV + skv = f{\v\^)v with /(O) = 0. 
In excitable media one has Ref(\v\'^) < 0, which provides the saturation of the instability. On the other hand, in 
CGLe the sign of Ref{\v\'^) is opposite and destabilizes the core according to simulations. Thus one has an alternative 
scenario of the meandering instability. The scenario appears to be generic and is not destroyed by small perturbations 
of the CGLe. 

Now going to larger e one finds that Rck increases with e and finally changes sign at a value £c. The result obtained 
from extensive numerical simulations is shown in Fig. 13. 

The acceleration instability can also be interpreted as the destabilization of a localized core mode similar to the 
situation of the standing hole solution in the Id CGLe (Popp et al, 1993, Sasa and Iwamoto, 1993; Chate and 
Manneville, 1992, see Sec. III). 

The case of negative diffusion ( £ < 0, |£| <C 1 ) is important especially for lasers with negative detuning (Oppo, 
1991; Staliunas, 1993; Lega et al, 1994). Then higher-order corrections to the diffusion (fourth derivative) have to be 
included to stabilize the short-wave instability. The resulting complex Swift-Hohenberg equation is 

dtA = A + lAA - e{ql + AfA - (1 + ic)\A\^A. (49) 

This equation exhibits in some range of parameters £, qc a similar core instability. In a certain parameter range one 
has stable meandering of the spiral (Aranson and Tsimring, 1995; Aranson, Hochheiser, and Moloney, 1997) 

C. Dynamics of vortices in the GLe, NLSe, and for 6 = c 

1. Dynamics of vortices in GLe 

In the GLe the spirals become vortices, i.e in Eq. (45) one has a; = and ip = Q. F{r) is a monotonic function and 
F ^ ar for r 0, F^(r) 1 — 1/r^ for r — > oo. 

A more general isolated moving vortex solution with wavevector Q away from band center is described by : 

A = B{r- vt) exp[iQr] (50) 
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where B fulfills the boundary conditions |_Bp 1 — for r — > oo with the appropriate phase jump of ±27r. ^ 

For large systems such that Rv 1, where R is the system size or the distance to another defect, the drift velocity 
V is perpendicular to the background wavenumber Q, such that the free energy decreases. The speed v is given by 
the expression 



with the constant vq = 3.29 (see Bodenschatz et al, 1988b, 1991, Kramer et al, 1990). Here U is the 7r/2 rotation 
matrix 



The term in the brackets on the rhs has been fitted to numerical results (sec Bodenschatz et al. 1991). The formula can 
be used throughout the stable range and is consistent with experiments of Rasenat et al (1990) in electroconvection 
in nematics. In small systems vR <C 1, i? » 1 the velocity is given hy v = 2Q/ log(_R/^o) with = 1-13. 

The above relations were rederived by Pismcn and Rodriguez (1991), Rodriguez and Pismen (1991), Ryskin and 
Kremenetsky (1991), and in part by Neu (1990a) using matching asymptotics techniques. 

The result shows that for an isolated defect the limits of the system size going to infinity and wavenumber mismatch 
(or defect velocity) going to zero are not interchangeable. This is analogous to the 2D-analog of Stokes' law (drag of 
an infinite cylinder through an incompressible fluid), where one also has a logarithmic divergence of the mobility at 
vanishing velocities. The similarity was pointed out by Ryskin and Kremenetsky (1991). It is also useful to note that 
the force exerted on dislocations in a phase gradient expressed by a wave vector Q is the analog of the well-known 
Peach-Koehler force on dislocations in crystals under stress. 

In the following we consider linearly stable periodic states, where < Q\. All periodic states, except Q = 
(band c{nitor), are metastable. Evolution to the band center can occur by the motion of vortices, once they have been 
nucleated by finite fluctuations or disturbances (see Bodenschatz et al, 1988b, and Hernandez-Garcia et al, 1993, for 
a discussion of the problem of homogeneous nucleation of vortices). 

The asymptotic interaction between two vortices can be described approximately by Eq. (51) with the r.h.s. replaced 
by the gradient of the resulting phase of other vortex at the center of the vortex. Since the gradient is perpendicular 
to the line connecting the vortex cores, the resulting force is directed along this line. As a result, oppositely-charged 
vortex attract each other and eventually annihilate, whereas like-charged repel each other. 

The analysis of an isolated moving vortex shows that the deformation of the roll pattern decays in front exponentially 
over a distance R ~ 1/v, while behind the dislocation the decay is proportional to R~^/^. Thus, if the background 
wavenumber is nonzero and two vortices with opposite topological charge arc approaching each other, the velocity 
will be constant for R < 2/v ~ 2/Q. Subsequently the motion will accelerate and eventually the attraction will 
dominate over the Peach-Koehler force and lead to annihilation. ^ In experiments by Brawn and Steinberg (1991), 
the annihilation process was studied in detail near the band center (Q <C 1) where accurate determination of Q is not 
possible. The quantitative comparison with theory is shown in Fig. 14 where the distance L between two vortices of 
opposite polarity approaching each other along a straight line versus time T is shown for different Q in scaled units, 
(Bodenschatz et al, 1991). At time T = the vortices annihilate. The solid line gives the numerical results, and the 
squares, circles and diamonds are the experimental data. 

The case Q = 0, where the Peach-Koehler force vanishes so that t; — > for L — > oo, deserves special attention. The 
analysis of Bodenschatz et al (1988)b, leads to vL = 2/C, where C ^ ln(L/2.26) for <C 1. Clearly for this to 
hold L needs to be exponentially large. Otherwise one is caught in the intermediate region tiL ~ 1 where no accurate 
relation for C was found. For that reason the so-called self- consistent approximation was proposed by Pismen and 
Rodriguez (1990), Rodriguez and Pismen (1991), Pismen (1999): 



®lt is important to mention that such a problem is not well defined for spiral waves (6 — c ^ 0). Since the spiral waves are 
active sources they can in general not coexist with the plane waves. Depending on the relation between the two frequencies 
either the spiral would invade the entire plane or the plane wave would push the spiral away. For details sec Sec. IV. F 

®In this discussion we have tacitly considered the situation where the Peach-Koehler force and the force due to interaction 
are parallel to each other. If this is not the case the vortices will move on curved trajectories. 




(51) 




ln{vo/v) = exp{±vr/2)[Ko{vr/2) ± Ki{vr/2)],vo = 3.29 



(52) 
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where the + sign corresponds to a hke-charged pair, the — sign to an oppositely-charged one, and Kq and Ki are 
modified Bessel functions. The formula is not very accurate, as shown by comparison with full simulations (Weber, 
1991, 1992a). 

Allowing for a very small wavevector displacement Q, one finds v{L + 1/Q) = 2/C. One can then get to the limit 
vL:^ 1 and then, for sufiiciently large systems, C = \n{3.29/v). 

For small distances, shortly before annihilation of an oppositely charged pair, the gradient terms in the GLe become 
dominant and one may argue that the dynamics is governed by self-similar solutions of the diffusion equation so that 
L ~ a/T, which appears consistent with the numeric results. 



2. Dynamics of vortices in the NLSe 

Although the stationary vortex solution of the NLSe coincides with that of GLe, the dynamic behavior of NLSe 
is very different from that of the GLe and is similar in the limit of small velocities, largo separation, and properly 
quantized circulation to the motion of point vortices in an ideal incompressible fluid (Lamb, 1932). The problem 
of the vortex motion in the context of the NLSe was first considered by Fetter (1966), and then many researchers 
rederived this result more accuracy (sec, e.g. Creswick and Morrison, 1980; Neu, 1990b; Lund, 1991, Pismen and 
Rubinstein, 1991a, Rubinstein and Pismen, 1994). 

Introducing the amplitude-phase representation A = i?exp[i^] and the "superfluid density" p = as well as the 
superfluid velocity V = —2Vd one obtains from Eq. (6) an Euler equation 

atV+(V-V)V = VP (53) 

and the continuity equation for the "density" of superfluid 

^tp-l-divpV = 0. (54) 

Here P = 2(— V^y^/y^ -|- p — 1) is the effective "pressure", whereas the first term is called "quantum pressure". 

The quantum pressure is irrelevant for large-scale perturbations of the density, i.e. V^p ^ p. In the Euler equation 
incompressibility approximation p = po = const ^ is reliable for the velocities much smaller than the sound velocity 
Cg = \/2po- One expects the analogy with classical incompressible vortex tubes in a liquid to hold for the well- 
separated vortices moving with the velocity u <C c^. Thus, according to classical hydrodynamics ( Lamb, 1932), one 
obtains that the oppositely charged vortex pair drifts with the velocity v = 2/L normal to the line connecting the 
cores of the vortices, and like-charged pair rotates around the center of the symmetry, where L is the distance between 
vortex center. This result up to higher order corrections holds for well separated vortices {L ^ 1) in the NLSe, as it 
was rederived later on by Neu (1990)b, Pismen and Rubinstein (1991), and by Rubinstein and Pismen (1994), using 
matching asymptotic technique. Numerical investigation of vortex dynamics in the NLSe was carried out also by Nore 
et al (1993), Frisch et al (1992), Abraham et al (1995), Josserand and Pomeau (1995). 
For A'' well-separated vortices located at the points Rj = {Xj,Yj) one has the more general formula 

^ 2yn,^^^-^^ (55) 

According to Eq. (55), obtained from the corresponding incompressible Euler equation, the vortex velocities diverge 
as the intervortex distances vanish. However, the above expression becomes incorrect when the inter-vortex distance 
d becomes of the order of the coherence length d ~ 0{1) (one in our scaling). For small separations the interaction of 
two like-charged vortices can be considered as a small perturbation of a double-charged vortex. Aranson and Steinberg 
(1995) have shown that for d ^ the velocity of vortices vanishes, and the frequency of rotation Q for like-charged 
vortex pair approaches a constant as d ^ (the "classical formula" (55) results in infinite frequency). Similarly, for 
an oppositely charged pair the velocity remains finite for d ^ Q. 

An important problem in the context of NLSe is the nucleation (or generation) of vortices from zero-vorticity flow. 
Frisch, Pomeau, and Rica (1992) considered the stability of flow passing an obstacle and found nucleation of vortex 
pairs above some critical velocity. Yet another mechanism of vortex generation is related to the instability of quasi- ID 
dark solitons, predicted by Kuznetsov and Turitsyn (1988). Josserand and Pomeau (1995) studied the instability of 
dark solitons numerically and found that the nonlinear stage of this instability results in creation of vortex pairs. 



23 



3. Dynamics of vortices for 6 = c 



The general case b = c, including the NLSe limit, can be related to the GLe limit by generalizing the ansatz (50) to 

A = B{r-vt)exp[iQr-ibt]. (56) 

Eq. (1) can be written as 



1 



62 



vV-B = AB + 2i{Q 



2(1 + 62) 



V VB+(1-Q^)B-|S|^B. 



(57) 



After neglecting the higher-order term Q^B, one observes that Eq. (57) coincides with the corresponding equation 
for 6 = if one replaces Q by Q + tttttt^v 



2(1+6^) 

As a result we may extend the expression (51) for arbitrary b leading to 



^Vlog((l+62)^;o/^;) =U 



(l + 62)Q + -v 



(58) 



One obtains oblique motion with respect to Q with the angle depending on b and Q. In the NLSe limit the motion 

becomes parallel to Q. Concerning the interaction between vortices one obtains oblique repulsion or attraction: for 
6^0 two vortices never move along the line connecting their centers but instead spiral with respect each other 
(Pismen, 1999). 



D. Dynamics of spiral waves for b ^ c 

1. General 

As we mentioned above, in the generic case of 6 ^ c the asymptotic interaction of topological defects (spirals) is very 
different from the corresponding dynamics of vortices in GLe and in NLSe. The problem of spiral wave interaction 
was considered by Biktashev (1989), Rica and Tirapegui (1990), (1991)a,b, (1992), Elphick and Meron (1991), Aranson 
et al (1991),(1993)a, Pismen and Nepomnyashchii (1991)a,b. The approach of Rica and Tirapegui (1990),(1991)a and 
Elphick and Meron (1991), explored a type of global phase approximation which in essence gave long-range interaction 
^ 1/r between the spirals, in analogy to the cases treated above. The results were in obvious disagreement with 
numerical simulations done by Aranson et al (1991), revealing exponentially weak interaction between the spirals. 

A more adequate solution of the problem was obtained later by Aranson et al (1991), and Pismen and Nepom- 
nyashchii (1991)a,b, on the basis of the phase diffusion equation (sec also Biktashev, 1989). The results demonstrated 
an exponentially weak asymptotic interaction, although they failed to describe the numerically observed bound states 
of spiral waves. The abovementioned approach is adequate for the case \b — c\ <^ 1. Biktashev (1989) considered a the 
very restrictive problem of spiral motion in an almost-circular domain. Due to some unphysical boundary conditions 
the results of Pismen and Nepomnyashchii (1991), (1992), were at conflict with numerical simulations (Aranson et al 
1993a). 

In later papers Rica and Tirapegui (1991)b,(1992) attempted to improve their interaction approach by combining it 
with that of Aranson et al (1991). However, in our view, the results remain unsatisfactory because of an uncontrolled 
perturbation technique. 

A consistent theory with quantitative predictions of asymptotic spiral wave interaction including bound-state for- 
mation was presented by Aranson et al (1993)a,b using the matching perturbation method. The idea of the method 
is straightforward: in the full solutions for a spiral pair (or a spiral and a wall, or more complicated aggregates of 
spirals) the spirals move with certain velocities, and thus solutions exist only with the "correct" velocities. Such 
solutions may be constructed approximately by starting with isolated spirals, each one restricted to the region in 
space filled by its emitted waves and moving with (small) velocities to be determined. For a symmetric spiral pair 
one simply has two half planes. In a first step the corrections are assiimed to be determined to sufficient accuracy 
by the linearized (perturbational) problem with boundary conditions that take into account the neighboring spiral 



The calculations can be generalized to arbitrary Q. 
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(or wall etc.). The veloeity eomcs in as an inhomogeneity. The occurrence of bound states can already be seen from 
the behavior of stationary perturbations of the asymptotic plane waves emitted by the spirals, as given by Eq. (25). 
Bounds states occur in the oscillatory regime, outside the curves (c— b)/{l + bc) = ±Ccr, Ccr = 0.845 (one must choose 
the curves inside the Benjamin- Feir-Newell-stablc region). The oscillatory range is plotted in Fig. 12 (curve OR). 
Also included in Fig. 12 is the convective Eckhaus-stability boundary (curve EI, see Sec. II. D) and the boundary of 
absolute stability (curve AI, see Sec. II. E) for the waves emitted by the spirals. 

The inhomogeneities involve the velocity linearly and as a consequence the sohition for arbitrary (but not too small) 
distance can be expressed in terms of one inhomogeneous and a homogeneous solution. Once these are determined 
numerically and used to match the boundary conditions the velocity versus distance relation comes out. For b = 
(for simplicity) the resulting velocities of the spiral at distance X from a plane boundary (or, equivalently, a pair of 
oppositely-charged spirals at the points {±X, 0)) are of the form 



J /-0\/l~-Q^cxp(-pX) rri /n \ 

V. = Im{ ^^^^^ X n/Im{C./Cy), 

vy = ReiZR^^E^^^x-n - v.Re{C./Cy). (59) 

where Q is the spiral wavenumber, Vx is the component of the velocity along the line connecting spiral cores (radial 
velocity) and Vy is the velocity perpendicular to this line (tangential velocity) . The constants Cx and Cy are obtained 
numerically and the parameter ^ is derived from the linearized problem (Aranson et al, 1993a). The bound states 
correspond to the case oi Vx = 0. The bound state will drift with the velocity Vy. The equilibrium distance 2Xe 
between the spirals is found from the equation 

Since only p and /x are complex this gives 

Im[pXe + iJ,hiXe] = —<p + ttI (61) 

where Z = 1, 2, 3... and (j> = — axg[\/ {5 Cyp^^"^)]. Even values of I correspond to stable bound states, odd I to imstable 
ones. Usually the lowest bound state with Z = 2 is relevant. Approaching the boundary of the oscillatory range the 
equilibrium distance diverges. 

Generalization to arbitrary b can be done along the lines of the similarity transformation (16). The result is at 
leading order unchanged. 

Note that the symmetries A{x, y) = A{—x, —y) (like-charged) and A{x, y) = A{—x, y) (oppositely-charged) have 
nearly the same effect on the boimdary conditions at x = for large spiral separation. Thus for large separation 
X the interaction of like-charged spirals is similar to the interaction of oppositely-charged ones. The only difference 
is that for the like-charged case both components of the velocities of the spirals have opposite sign whereas for the 
oppositely-charged case the tangential components have the same sign. This causes the rotation of the spirals around 
the common center of the symmetry in the like-charged case. 

Like-charged spirals may form more complicated bound states or aggregates. In contrast to the two-spiral bound 
states, which are simply rotating with constant velocity, each spiral in the aggregate performs a more complicated 
motion (possibly nonperiodic) on the background of a steady-state rotation, perhaps similar to the "dance" of spiral 
aggregates observed in a variant of the Belousov-Zhabotinsky reaction (Steinbock and Miiller, 1993). Certain problems 
of interaction between spiral waves and hole solutions were considered by Bazhenov and Rabinovich (1993,1994). 
Interaction of spiral with external periodic perturbation and response functions of the spiral core were studied by 
Biktasheva et al (1998,1999,2000). 



'^'^Note that the additional term —ihwVA arising in the corresponding equations of motion for each spiral can be absorbed in 
the nonsingular part of the linear problem (Aranson and Pismen, 2000). 
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2. Comparison with results of numerical simulations 



Bound states of spirals are shown in Fig. 15. In Fig. 16a the dependence of the velocities on the spiral separation 2X 

is plotted for & = 0, c = 1 and compared with results from full numerical simulations. There is reasonable agreement, 
particularly for the radial velocity Vx- The first two zeros of Vx at 2X(. ft! 11.5 and 22.8 correspond to / = 1, 2 
in Eq. (61). Here the velocities are already extremely small. From the simulations no other bound state could be 
resolved. The equilibrium distance obtained from the theory for the stable boimd state (Z = 2 in Eq. (61)) is in a very 
good agreement with the results of the simulations of the full CGLe (see Fig. 16a). There is a discrepancy between 
the values of the tangential velocity Vy by a factor of about 1.5 (see Fig. 16b) although the functional dependencies 
on c are very similar. This discrepancy results from the fact that the shock is not described well in linear treatment. 

More accurate estimates can be obtained by treating the shock fully nonlinearly. Far away from the core the motion 
of the spiral can be neglected and one may consider the stationary shock produced by the waves emitted by the spirals. 
Thus one has to solve the stationary CGLe with boundary conditions A ^ r cxp[z6'] for a; — > X, and dxA = at x ^ 0. 
This is a rather complicated 2D problem. Numerically one could extract more accurate values of the constants C\^2m 
characterizing reflection at the shock, see Fig. 16b. 

3. Interaction in the monotonic range. 

In the case < c < Ccr there exist two real positive roots of Eq. (25). Using the analogous numerical procedure one 
can also determine the constants Cx.y, but for small c it is technically very difficult in this form. The results can be 
simplified considerably for the case c ^ and \cQ\X 1. Then one can neglect the coefficients C2„ because for c ^ 
one has < pi w -2cQ < l,p2 « v^,/ii and -1/{2Q'^). One obtains Vx,y « C^^ exp(-2|c<3|X)/v^, 

which exhibits explicitly the exponential decay of the interaction and reproduces the earlier analysis using a phase 
diffusion equation (Aranson et al, 1991; Pismen and Nepomnyashchii, 1991,1992; Biktashev, 1989). 

Numerical simulations for oppositely charged spirals for 6 = and c = 0.5 indicate asymptotic repulsion (Aranson 
et al, 1993a). This result can bo inferred already from the work of Biktashev (1989), who used a phase-diffusion 
equation and asymptotic matching to treat the interaction of spirals with a boundary. It appears to be at conflict 
with the works of Pismen and Nepomnyashchii (1991)a,b where in the limit c — > matching with the internal solution 
was done analytically and asymptotic attraction was found. The repulsive range is expected to move to larger X 
roughly as |cQ|~^ for smaller c. In this way the crossover to the long-range attraction within the GLe is recovered. 
Interaction becomes attractive at smaller distance leading to final annihilation of the spiral pair. 

Like-charged spirals for c < Ccr have repulsion at large distance (as in the oppositely-charged case) and also at small 
distance. So it is quite clear that the interaction is repulsive everywhere. 

E. Interaction of spirals with a inhomogeneity 

Spirals may interact with inhomogeneities of the medium. This problem is particularly interesting in the context of 
Bosc-Einstein condensation as described by the NLSe in a parabolic potential (for a review see Dalfovo et al, 1999). 
Weak inhomogeneities are sometimes included into the CGLe in the form 



with an appropriate function u{r) . An interesting problem is the drift of a spiral in the gradient created by a localized, 
radially symmetric inhomogeneity (see Staliunas, 1992, Gil et al, 1992). The spiral may be trapped (pinned) by the 
inhomogeneity or perform a stationary rotation at some distance from the inhomogeneity as observed in optical 
systems (Arecchi, 1990, 1991; Brambilla et al, 1991). 

Using the method described above one arrives at the following expression for the spiral velocity due to interaction 
with weak axisymmetric inhomogeneity (Aranson et al, 1993c) 



where r is the distance from the inhomogeneity to the spiral core, v is the spiral core velocity and the function G 

is obtained by a convolution of the inhomogeneity profile ^{r) with the corresponding solution of linearized CGLe. 
Equation (63) is the analog to the motion of a " massless particle" in a radially-symmetric field. Numerical simulations 



dtA = {1 + i^{r))A + {1 + ib)AA -{1 + ic)\A\^A. 



(62) 



V = G(r) 



(63) 
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with CGLe shows that Eq. (63) describes the motion of the spiral fairly well for not too large values of see Fig. 
17a. 

For large |6| the complex trajectories of the spiral core are not captures by Eq. (63), see Figure 17b. In order to 
describe the numerical simulations one has to include the effective "mass" term into the equation of motion: 

f+.^v = G (64) 

where K is the friction tensor (or mobility tensor), and e = 1/6, see Sec. IV. B. 2. Although these features were 
obtained in the stable range Rck > {e > Ec), the external forces created by the inhomogeneity can excite the weakly 
damped core mode. This can explain the complex motion of the spiral core observed in the presence of obstacles 
(Sepulchre and Babloyantz, 1993). 

Spiral motion in a slowly varying (on the scale of the coherence length) inhomogeneity was considered by Hedrey 
et al, (1999). In this situation additional simplifications are possible. A related problem is the motion of the spiral 
in the presence of small thermal noise. As was shown by Aranson et al (1998), the spiral core has finite mobility and 
performs Brownian motion. 



F. Symmetry breaking 

Symmetric bound states of spirals are not necessary stable. There is numerical evidence that symmetric bound 
states, after a sufficiently long evolution, spontaneously break the symmetry such that one spiral begins to dominate, 
pushing away other spirals (Lcga, 1990; Weber et al, 1991a, b; Aranson et al, 1993b). Synmictry-brokcn states arc 
often produced directly from random initial conditions (Aranson et al, 1993b) . To understand the symmetry- breaking 
instability one must consider the perturbation of the relative phase, or the frequency ui of the waves emitted by each 
spiral, caused by the interaction with the other spiral. Indeed, from the analysis of CGLe it is known, that the shock 
(or sink) where two waves with different frequencies uii collide moves in the direction of smaller frequency, due to 
conservation of phase. This means that after sufficiently long time only the larger frequency (or equivalently the 
larger wavenumber \Q\, because of the dispersion relation u — c(l — Q^) + 6(5^) dominates in a bounded system. The 
velocity of the motion is given by ~ (c — 6)(Qi + Q2) = ^{vgi + Vg2) where Vg = duj/dQ is the group velocity of 
a plane-wave state. If due to the interaction the frequencies of the spirals become different, one can expect a drastic 
breaking of the symmetry of the system. This effect appears to be very important for the generic long-time evolution 
of large systems containing spirals (vortex glass). 

The symmetry breaking instability of spiral pairs can be understood from the dynamics of the relative phase of the 
spirals 4> = <t)i — (j)2- Aranson et al (1993b) derived an equation governing the phase difference between two spirals 
separated by the distance 2X 

= (^^"^^^^^^^""^0 ^^'^^(i)) /^-(^-^o ) - (65) 

where Cio,Co are the constants obtained from the linearized CGLe. The last term represents the lowest-order term 
of an expansion in (f). The constant C determined numerically from Eq. (65) turns out to be positive for the (first) 
symmetric bound state (sec Fig. 18), so that the state is unstable with respect to 0. According to Aranson et al 
(1993b), bound states are unstable with respect to symmetry breaking in the whole oscillatory range. Due to the 
fact that C, is very small, the symmetric bound state is rather long lived. As a result of the symmetry breaking only 
one "free" spiral will remain, whereas the other spiral is pushed away to the boundary. Depending on the boundary 
conditions the second spiral will finally either annihilate at the boundary (non-flux boundary conditions, i.e. zero 
normal derivative on the boundary), or, with periodic boundary conditions the defect will persist for topological 
reasons, but reduced to a sink and enslaved in the corner of the shock structure of the free spiral (edge vortex) (see 
Fig. 19). This leads to an asymmetric lattice of topological defects, which appears to be stable in the oscillatory case. 

The absence of symmetry breaking in the limit |c — 6| ^0 and |(c — b)Q\X ^ 1 can be extracted already from 
the paper of Hagan (1982). Prom this work one finds that decreasing the radius increases the spiral wavenumber, 
which ultimately means stability of the shocks. The stability of symmetric states together with the repulsion of the 
spirals indicates the possibility for the existence of symmetric ( "antiferromagnetic" ) lattices of spirals (or "Wigner 
crystal"), i.e. lattices made up of developed spirals with alternating topological charge. Such lattices were obtained 
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in numerical simulations in this range by Aranson et al (1993b). They appear to exist in the whole non-oscillatory 
range 

Since for large separation X the properties of like-charged spirals are analogous to those of oppositely-charged 
spirals, one expects the same mechanism of symmetry breaking. Indeed for |6 — c| above the critical value such a 
breaking was observed in numerical simulations. 

Certain aspects of the spiral pair dynamics were studied numerically by Komineas et al (2001). It was shown that 
for c > Ccr (symmetry breaking range) the spirals interchange partners and form new pairs. For c < Ccr (monotonic 
range) symmetric spiral states were found. 

G. Vortex glass 

A fascinating problem concerns the effect of the symmetry breaking on the long-time behavior of large systems that 
evolve from random initial conditions in the parameter range where one has the oscillatory interaction and spirals. 
Numerical simulations, carried out by Huber et al (1992), Chatc and Manncvillc (1996), have shown a remarkable 
phenomenon: from an initially strongly turbulent state eventually spirals with a developed far field evolve. The spirals 
grow until the entire space is filled (see Figure 20), each spiral occupying a certain domain. The domains have various 
sizes and typically have a four or five-sided near polygon structure. Locally each domain boundary, represented by 
a shock, is nearly hyperbolic. The shocks often contain "enslaved", or edge vortices. This state appears to persist 
indefinitely, although sometimes domain boundaries break and edge vortices annihilate. The state was called "vortex 
glass" by Huber et al (1992). Some aspects of relaxation to the vortex glass were considered by Braun and Feudel 
(1996). 

The existence of the "vortex glass" is connected with the nonmonotonic dependence of the spiral frequency on the 
domain size. As one sees from Figure 21 the spiral frequency oscillates as the radius of the domain increases and 
approaches the asymptotic value for r ^ oo. (The computations were done for a radially symmetric domain, but 
we expect the result to be qualitatively correct for domains with arbitrary convex shape.) As a result, one has sets 
of domain radii r,; corresponding to the same frequency. Neighboring domains can only coexist without (further) 
"symmetry breaking" if they correspond to the same spiral frequency. Ultimately one can conceive of a network of 
fully synchronized domains of various shapes and sizes. These are likely the ingredients necessary to explains the 
variability of the domain sizes observed in the quasi-stationary vortex glass. The vortex glass is presumably the global 
attractor also in the full convectively stable oscillatory range. 

Bohr et al (1996), (1997) considered the shapes of spiral domains in the vortex glass. They found that they differ 
from Voronoi Polygons and obtained the form of the domain boundaries from the condition of phase continuity across 
the shocks. In particular, far away from the center of an unperturbed spiral, the phase is given by 4>i = ±9i + Qri + Ci, 
where rj, 6i are the polar coordinates measured from the center of the spiral and C, is the phase constant of the spiral. 
If one assiimes that the distance between two neighboring spirals is much larger than the wavelength 211 /Q, one 
obtains simply a hyperbola describing the shock shape: — rj = {Cj — Ci)/Q. This simple formula reproduces the 
structure of the spiral domains with high accuracy (see Figure 22). 

H. Phase and defect turbulence in two dimensions 

On a rough scale, two types of turbulent behaviors can be identified in 2D: phase and defect turbulence. Defect 
turbulence is characterized by persistent creation and annihilation of point defects. By contrast, in phase turbulence 
no defects occur. In the system with periodic boimdary conditions the total phase gradient across the system (the 
"winding rnimber") is conserved. The simulations of Manneville and Chate (1996) provide some evidence that phase 
turbulence breaks down in the infinite size, infinite time limit. However, it could be that the transition is characterized 
by a "sharp" bifurcation scenario similar to the situation in ID (see III.D.l). 



^■^ Aranson et al (1993b) have found stable square lattices in rather small system and for relatively small number of spiral 
(4x4). There is no guarantee that similar lattices are stable in arbitrary large systems. 
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1. Transition lines 



The most detailed survey of various regimes occurring in 2D is given by Chate and Manneville (1996) and Manneville 
and Chatc (1996). According to Chatc and Manneville (1996), the transition from vortex glass to defect turbulence 
starting from random initial conditions occurs at the numerically determined line "T" , see Fig. 12. The transition 
occurs somewhat prior to the absolute instability limit given from the linear stability analysis of plane waves emitted 
by spirals. However, starting from carefully prepared initial conditions in the form of large spirals one can approach 
the absolute instability limit. 

Before the line "T"one finds transient defect turbulence which finally exhibits spontaneous nucleation of spirals 
from the "turbulent sea". At the line "T" the nucleation time presumably diverges. Before the line "T" the entire 
space will finally be filled by "large spirals" separated by shocks, forming a vortex glass state. 

As in ID, persistent phase turbulence exists also in 2D between the Benjamin- Feir-Newell line and the line "L" 
(Chate and Manneville, 1996; Manneville and Chatc, 1996). The range is somewhat smaller than in ID. Beyond 
the line "L" defects are created spontaneously leading to defect chaos. Actually, in 2D, phase turbulence is always 
metastable with respect to defect turbulence or vortex glass. 

2. Spiral break-up 

Spiral break-up was observed in experiments on chemical oscillatory media by Ouyang and Flesselles (1996), Zhou 
and Ouyang (2000), in simulations of CGLc by Chatc and Manneville (1996) and model reaction-diffusion systems 
by Bar and Or-Guil (1999). The instability likely originates from the outer region of the spiral and thus is expected 
to occur beyond the line "T". Quenching a large spiral to this domain will usually trigger a breakup scenario. On 
the other hand, a careful adiabatic procedure allows, in principle, to avoid breakup and to preserve one large spiral 
up to the (linear) absolute instability threshold. (In the presence of noise, the convective instability actually imposes 
a maximal noise-dependent radius iinoise-) Not too far from the absolute instability limit one is left with a smaller 
spiral surrounded by strong chaos (Fig. 20, left) whose well-defined radius R does not depend on system size and/or 
distance to the boundaries but vanishes as one approaches the absolute instability threshold. 

In the experiments and in some model systems, a spiral need not be introduced initially, as the chaotic phase is only 
metastable to the spontaneous nucleation of spirals whose radius grows up to R. Thus, the asymptotic configuration 
on a long time scale is actually a quasi- frozen cellular state (Fig. 20, right). 

The above concept of spiral break-up in large systems is in contrast to arguments by Tobias and Knobloch (1998), 
who stated that the spiral wave breakup occurs in the regime of absolute Eckhaus instability via "a globally unstable 
wall-mode confined to the outer boundary" whose front structure is at the origin of the stable "laminar" spirals 
immersed in a turbulent sea. Whereas a stationary front can indeed be associated with the absolute stability limit 
(see Sec. II. E) it can only give an upper limit to the existence of the convectivcly unstable state. 

In the break-up scenario studied by Zhou and Ouyang (2000) spirals emitting modulated waves in the presence 
of stable meandering of the core were observed. This behavior can be understood in terms of spatial amplification 
of periodic perturbations due to core meandering in the regime of convective instability of the background waves 
(Brusch et al., 2001). Although the saturated meandering does not occur in the framework of the CGLe, a properly 
perturbed CGLe may exhibit both, saturated meandering and convective instability, see Aranson, Hochheiser and 
Moloney (1997). 

3. Defects statistics 

Defect turbulence is the most chaotic state in 2D. It is characterized by exponential decay of correlations, with 
short correlation lengths and times. The density of defects varies with b and c. However, it can be argued that, at 
least in the most studied region with b, c of order one, the defects do not play the role of "particle-like" excitations in 
this case. Indeed, the defects in the defect turbulence are very different from spiral waves since they don't emit waves. 
They behave as passive objects and are merely advected by the surrounding chaotic fluctuations. According to Chate 
and Manneville (1996) "amplitude turbulence" is a more appropriate name for such spatio-temporally chaotic state. 

One may argue that the stationary distribution for the number of defects is described by p{n) ~ exp[— (n — n)^/fi] 
where n is the average number of defects. The formula follows in the limit n ^ oo from treating defect pairs as 
statistical independent entities (Gil et al, 1990). Note that this assumption is strictly valid only for large domains. 
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Otherwise the statistics is influenced by defects entering/leaving the subsystem. When these processes dominate the 
exponent acquires a factor i. 

Egolf (1998) extracted the degrees of freedom associated with defects from those of phase fluctuations by using 
the concept of finite-time Lyapunov dimension. He found that each defect "carries" from one to two degrees of 
freedom. Although one may argue that the number of defects can be a convenient characterization of some types of 
spatio-temporal chaos, the method of separation does not appear fully convincing because the relation between the 
finite-time dimension and number of defects is not examined for different system sizes and duration intervals. 

It is interesting to mention that a similar analysis performed by Strain and Greenside (1998) for a reaction diffusion 
system results in a much higher dimension per/defect, namely between 3 and 7. 

Mazcnko (2001) studied defect statistics in the CGLe in the "defect-coarsening regime" (presumably the region of 
small b and c and mean distance between defects smaller than the screening length 1/ {\b—c\Q)) using an expression for 
the defect velocity similar to that of Rica and Tirapegui (1990), (1991a), which is not applicable for defects separated 
by more than the core size of 0(1). 

4. Core instability and spiral turbulence for large 6 

One expects that in the full core- unstable range s = 1/b < Sc, a state with persistent defects should be typically 
spatio-temporally chaotic. 

In the monotonic range (Fig. 13, above curve OR), this turbulence is characterized by fast motion of the defects 
and collisions, whic;h often do not result in annihilation, in contrast to the usual defect chaos discussed above. In Fig. 
23, the number n of defects as a function of time is shown for a fairly large system (150 x 150). Apart from the rapid 
fluctuations due to creation and annihilation there is an extremely slow decrease of n. Crossing Sc in the monotonic 
range the disordered state appears to persist (or is at least very long lived). This indicates a hysteretic behavior, 
which is to be expected from the suberitieal character of the core instability. 

In the oscillatory Eckhaus stable range (Fig. 13, below curve OR) the behavior is drastically different. Starting 
from random initial conditions one first has the evolution towards a vortex glass as in the range e > £<;, see Fig. 
23 and Sec. IV. G (Aranson et al 1993b). However the spirals are unstable with respect to acceleration resulting in 
continued dynamics of the spiral and shocks. Additional spirals are created very rarely. One might call this a "hot 
vortex glass" . 

When the Eckhaus instability sets in at the curve EI of Fig. 13 the perturbations produced by the accelerating 
core of the dominant spiral are amplified away from the core due to the convective character of the instability. When 
some critical level is exceeded, the state looses stability and many new defects are created throughout the cell. Then 
the process repeats (Figs. 23 and 24). Such phenomena are very similar to spatio-temporal intermittency of holes 
observed in the ID CGLe (Chate et al, 1994; Popp et al, 1993, see Subsec. III.D.4) and might be called defect-mediated 
intermittency. In a large cell one expects to have such processes developing independently in different places of the 
cell, so one has persistent chaotic bursts (or spots) on the background of growing spirals. 

Below the curve ST in Fig. 13, which presumably is the continuation of the T curve in the range s > Sc (Chate 
and Manneville, 1996), the "strong" chaotic characteristic of the intermittent bursts becomes persistent. The state 
is similar to the usual defect chaos. The curve ST lies somewhat above the limit of absolute instability, curve AI(see 
Sec. II. E). For e — > the Al curve tends to c « —1.2. The curve was determined by simulating Eq. (46) with the 
restricted class of functions (45) and boimdary conditions ^4(0) = Q,drA{L) = 0,L ^ 1, which is effectively a ID 
problem. The curve ST determined in this manner is consistent with 2D simulations (Aranson et al, 1994). 

V. DYNAMICS IN 3D 
A. Introduction 

The 3D analog of the 2D vortex or spiral wave is called vortex filament or scroll wave. The point singularity of the 
phase of the complex function A at the center of the spiral becomes a line singularity in 3D. The filaments can be 
open (scrolls), closed (vortex loops and rings), knotted or even interlinked, twisted, or entangled. Depending on the 
parameters of the CGLe the scroll wave can be stable or develop some instabilities. Remarkably, 3D vortices can be 
highly unstable even in the range of parameters where their 2D analog is completely stable. 

CGLe has a stationary solution in the form of a straight vortex with twist 



30 



A{r, e, z) = F{r) exp i[a;t ±6 + ip{r) + k^z] , 



(66) 



Here the axial wavenumber kz characterizes the twist. Curved vortex hnes are non-stationary. In most of the cases 
vortices untwist, and the solution with fc^ = is the most stable one. 

Scroll waves have been observed experimentally in slime mold (Siegert and Weijer, 1991), heart tissues (Gray 
and Jalife, 1996), gel- immobilized Belousov-Zhabotinskii (BZ) reaction (Vilson et al, 1997). Long-lived entangled 
vortex patterns in three-dimensional BZ reactions were observed by the group of Winfree using optical tomography 
techniques, (Winfree et al, 1996). Complex vortex configurations have also been observed in numerical simulations 
of reaction-diffusion equations (Winfree, 1995, Biktashev, 1998, Fenton and Karma, 1998a,b; Aranson and Mitkov, 
1998; Qu, Xie and Garfinkel, 1999). 

Theoretical investigation of scroll vortices in reaction-diffusion systems was pioneered by Keener and Tyson, 1990, 
1991, who derived the equation of motion for the filament axis. In particular, it was found that vortex rings typically 
shrink with a rate proportional to the local curvature of the filament, leading to collapse in finite time. The existence 
of non- vanishing vortex configurations and expansion of vortex loops, observed also in numerical simulations of 
reaction-diffusion equations, was interpreted as "negative line tension" of the vortex filament (Biktashev, 1998). 

In order to characterize the motion of vortex lines in three-dimensional space let us consider a curve C at any 
moment of time i in a parametric form X(s, t), where s is the arclength. At any point of the curve a local orthogonal 
coordinate basis can be defined {Frenet trihedron), yielding the Prenet-Sorret equations (see e.g. Pismen (1999)) 

Xg = 1, Ig = nn, Ug = —kI + rb, = — rn (67) 

where 1, n, b are tangent, normal and bi-normal unit vectors, and k and r are curvature and torsion of the curve, 
correspondingly. 



B. Vortex line motion in NLSe 

As in 2D, the analysis of vortex motion in the 3D NLSe is based on the analogy with the Euler equation for ideal 
fluids, (see Batchelor, 1967; Pismen, 1999). Following the analogy with the vortex lines in an ideal fluid, the local 
"superfluid velocity" v = (j) = argA can be found from the Biot-Savart integral 

F / R X dl 

where J> S/cj) = 27rF, F = ±1 is the vorticity of the line, R = x — X, R= |R|. The Biot-Savart integral describes the 
velocity far away from the vortex line but does not allow to compute the velocity of the lino since it diverges at the 
core. In order to get the velocity of line motion one need to perform matching of the vortex core field with the far 
field given by the Biot-Savart integral (Pismen and Rubinstein, 1991, Pismen, 1999). As a result of matching one 
obtains 

= Vs + FKb log — (69) 

«0 

where vj^ is the vortex drift velocity in the local normal plane, A is a constant dependent on the geometry of vortex, 
and oo ~ 1.856 for a single-charged vortex line. For a vortex ring of the radius R one has A = 8R. The last term 
expresses the drift along the bi-normal with the speed proportional to the curvature, the first term Vg accounts for a 
contribution from non-local induction of the Biot-Savart integral. Neglecting the first term one recovers the so-called 
localized induction approximation, often used in hydrodynamics ( for review see Ricca (1996)). The localized induction 
approximation is in fact a rather crude approximation for vortex motion in the NLSe. It has been extensively used 
due to its mathematical elegance. 

In the simplest case of the vortex ring of radius R one obtain from Eq. (69) that the ring drifts as a whole along 
the bi-normal, i.e. along the axis of symmetry, with the speed given by (Pismen, 1999) 

v=^log{R/Ro) (70) 

where Rq « 0.232 is the constant obtained from matching with the vortex core. Clearly, a similar result can be 
obtained for classical vortices in an ideal fluid with the difference that in the latter ones there is no a well-defined 
core. 
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C. Collapse of vortex rings In the CGLe 



In GLe one finds that a vortex ring (radius R) collapses. Indeed, substituting the ansatz 

A{r,0,z,t) = Ao{r- [ v{t')dt',z) (71) 
Jo 

{Aq is the 2D stationary vortex solution; r, 6, z are cylindrical coordinates) into the 3D GLe one derives 

-vdrAo = AqU - l^oH + d^Ao + dlAo + UrA^. (72) 

Replacing the explicit r dependence in the last term by the radius of the ring, one obtains from the consistency 
condition that 

dR 1 



Solving Eq. (73) one derives R(t) = \JRq — 2t, i.e. the ring collapses in finite time. Surprisingly, in the case of the 
GLe the analog of the localized induction approximation produces the correct answer. 

Gabbay et al (1997), (1998a) generalized this result for the CGLe where the ansatz Eq. (71) has to be generalized 
to include the curvature- induced shift of the filament wavenumber. They showed that the ring collapses in finite time 
according to the evolution law 

^ = -l±^. (74) 
dt R ^ ' 

In addition, there is no (at least, at first order in 1/R) overall drift of the vortex ring in the direction perpendicular 

to the collapse motion. The collapse rate (often associated with the "line tension") v ~ 1 + 6^ appears to be in 
reasonable agreement with simulations for not too large \b\. This corrects a previous erroneous estimate v = I + be 
(Prisch and Rica, 1992). 

For the evolution of the local twist of a straight vortex one obtains the Burgers equation (Gabbay, Ott, Guzdar 
1997; 1998b, Nam et al, 1998) 

dt<l> ={b- c)(l - k^){d,cj>,f + {l + be+{b- e)bkl)d,,(t> (75) 

where dzii = kz and fco is the asymptotic wavenumber of the 2D spiral solution. A more complicated equation was 
obtained for a twisted and curved vortex filament. 



D. Vortex nucleation and reconnection 

Vortex reconnection in NLSe was studied in relation to turbulence in supcrfluid liquid helium. Large-scale compu- 
tations were performed by Schwarz (1988) using the localized induction approximation. These computations give an 
impressive pictures of vortex tangles. However, it remains unclear if the particular features of the tangle are real or 
an artifact of the localized induction approximation. 

Obviously, vortex reconnection must be described by the full NLSe. Koplik and Levine (1993,1996) find in full 
numerical simulations of the NLSe that vortices reconnect when they approach within a few core lengths. Depending 
on conditions, the vortex rings may scatter, merge and break up subsequently, sec Fig. 25 

Gabbay et al (1998b) studied vortex reconnection in the CGLe. As a result of the interplay between two effects: 
motion of the filaments towards each-other due to attraction (two-dimensional effect) and opposite motion due to 
curvature (three-dimensional effect) a criterion for vortex reconnection was proposed. 

In the GLe vortex rings ultimately shrink. However, with an additional phase gradient j parallel to the ring axis 
(e.g. due to rotation of a superfluid or supercurrent in superconductors) the vortex rings may expand. The additional 
force is the analog of the Pcach-Kochler force on a 2D vortex in a background wavenumber Q = j. In this situation 
one has the following equation for the ring radius R: 

dR 1 



32 



where a k, 2j/\og{vo/2j) (sec Eq. (51)). Thus, if i? > l/{aj) the vortex ring will expand. 

Expanding vortex rings in the GLe were obtained in simulations as a result of nucleation after a rapid thermal 
quench by Aranson, Kopnin and Vinokur (1999). This problem was considered in the context of experiments in 
supcrfluid liquid ^He heated well above the transition temperature by absorption of neutrons by Rimtu et al (1996, 
1998). This experiment was designed to verify the fluctuation-dominated mechanism for the formation of topological 
defects in the early Universe suggested by Kibble (1976) and Zurek (1985) and elaborated in later work by Dziarmaga 
et al (1999), Antunes et al (1999) mostly in ID and 2D. 

Selected results are shown in Fig. 26. One sees (Fig. 26a-c) that without fluctuations the vortex rings nucleate 
upon passage of the thermal front. Not all of the rings survive: the small ones collapse and only the big ones grow. 
Although the vortex lines arc centered around the point of the quench, they exhibit a certain degree of entanglement. 
After a long transient period, most of the vortex rings reconnect and form an almost axisymmetric configuration. 

It turns out that fluctuations have a strong effect at early stages: the vortices nucleate not only at the normal- 
supcrfluid interface, but also in the bulk of the supercooled region (Fig. 26d-e). However, later on, small vortex rings 
in the interior collapse and only larger rings (primary vortices) survive and expand (Fig. 26f). 

E. Instability of weakly-curved filaments in the large b limit 

Aranson and Bishop (1997), Aranson let al (1998) have shown that the simple relation for the collapse rate Eq. 
(74) is violated in the large b limit, 6=l/e;^'l. As a result of an asymptotic expansion for e <^ 1 the equation of 
motion of the filament takes the form 

Otv + K[ev - nn] = 0. (77) 

where v is the velocity of the filament, and n the Ioc;aI curvature. The 2x2 matrix K corresponds to the "friction" 
matrix of the 2D spiral waves, see Sec. IV. B. 2. 

Note, that dropping the acceleration term in Eq. (77) one recovers the result of Eq. (74) for 6 ^ oo, since for the 
ring ujv = dtR, n = ~1/R. Restoring the original scaling r — s- r/Vb, one obtains dtR = —b^/R. However, since in 3D 
the local velocity in general varies along the vortex line, even small acceleration may cause severe instability, because 
the local curvature becomes very large. Moreover, deviation of the local velocity from the direction of the normal will 
lead to stretching and bending of the vortex line. Thus the acceleration term, which formally can be considered as a 
higher-order correction to the equation of motion, plays a pivoting (crucial) role in the dynamics of a vortex filament. 

1. Perturbation around a straight vortex 

An almost straight vortex parallel to the z axis can be parameterized by the position along the z coordinate: 

(Xo(z), lo(2))- Since in this limit the arclcngth s is close to z, the curvature correction to the velocity Kn is simply 
KW = (S^Xq, OIYq) = d-^r, where r = {Xq, Yq). Using Sjr = v, Eq. (77) reduces to a linear equation 

dtv + K[ev - d^r] = 0, (78) 

The solution can be written in the form r ~ exp[ikz + X{k)t], where A is the growth rate. We immediately obtain the 
following relation for A: 

A^ + x(eA + = 0. (79) 

where x = ^xx i iK^y is the "complex friction'' ( compare Sec. IV. B. 2). One may consider two cases: fc e and 
k ^ e. For A: <C e from Eq. (79) one obtains A = —ex + C>(fc^), i.e. recover the core instability of the 2D spiral. For 
k ^ e one derives A w ±^ —(Kxx ± iKxy)k. There always exists a root with a large positive real part: \ ^ k ^ e. 
Therefore, for finite k, the growth rate \{k) may significantly exceed the growth rate of the acceleration instability 
in 2D (corresponding to fc = 0): ReA — —eK^x- Hence, the "small-curvature" approximation considered above can 
be valid only for finite time. The fall off of the growth rate A at large k is not captured by the small-curvature 
approximation used here. 
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2. Numerical results 

Comparison of the theoretical results with numerical simulations is shown in Fig. 27. As an initial condition a 

straight vortex line with small periodic modulation along the z-axis was taken. As one sees from the figure, the 
growth rate indeed increases initially with k, and then falls off for large k. The theoretical expression (79) shows 
reasonable agreement with the simulations for small enough k. The growth rate at the optimal wavenumber exceeds 
the corresponding growth-rate of the acceleration instability (fc = 0) by more then two orders of magnitude. 

The long-time evolution of a perturbed straight vortex is shown in Fig. 28. As one sees from the figure, the length 
of the vortex line grows. The dynamics seems to be very rapidly varying in time, and the line intersects itself many 
times forming numerous vortex loops. The long-time dynamics shows, however, a saturation when a highly- entangled 
vortex state is developed and the total length of the line cannot grow further due to a repulsive interaction between 
closely packed line segments. The dependence of the line length on time is shown in Fig. 29. One can identify two 
distinct stages of the dynamics: first, fast growth of the length; second, oscillations of the line's length around some 
mean value. 

For small enough e, two distinct behaviors of the total vortex length depending on the value of c are observed. Above 
a critical value Cc corresponding approximately to the convective instability range of the 2D spiral (cc for e ^ 0), 
the total length approaches some equilibrium value and does not exhibit significant fluctuations. On the contrary, 
for c < Cc, the total length exhibits large non-decaying intermittent fluctuations around the mean value. Figure 30 
presents the snapshots illustrating the structure of the vortex field corresponding to the moments of maximum and 
minimum of the length. One sees, that in this situation some segments of vortex lines start to expand spontaneously, 
pushing away other vortex filaments and in such a way making substantial vortex-free holes around them. Then the 
instability takes over and destroys these almost-straight segments of filament, bringing the system back to a highly- 
chaotic state. This dynamics can be considered as a 3D spatio-temporal vortex intermittency, which is an extension 
of spiral intermittency discussed in Sec. IV. H. 4. For even smaller values, c < — 1, one has the transition to a highly 
chaotic state, which is an analog of "defect turbulence" in the 2D CGLe. In this regime small vortex loops nucleate 
and annihilate spontaneously. 

The evolution of a closed vortex loop is shown in Fig. 31. The simulations show that the 3D instability may prevent 
the ring from collapse, causing the stretching of the loop in the direction transversal to the collapse motion. However, 
small rings typically collapse, since then the instability described above docs not have time to develop substantial 
distortions of the ring. Even in this situation the ring exhibits a few oscillations of the radius. 

3. Limits of three-dimensional instability 

The previous analysis indicates instability of vortex lines in the limit e — > for all c. However, it cannot describe 
the boundary of the instability for increasing e. In order to obtain the stability limit one needs to perform a full 
linear stability analysis of a straight vortex solution, not limited to small k and e (Aranson, Kramer, Bishop, 1998). 
The linear stability analysis shows that the 3D instability persists substantially beyond the 2D core instability. The 
results are shown in Fig. 32. Moreover, the typical growth-rate in 3D is much higher than in 2D. 

F. Helices, twisted Vortices and supercoiling instability 

Close to the stability boundary of the 3D instability, the evolution of a straight vortex does not necessarily result in 
spatio-temporal chaos. In contrast, the simulations show that the instability saturates, resulting in a traveling helix 
solution (Aranson et al 1998) or superposition of two helices with opposite chirality (Rousseau et al 1998). Indeed, since 
the left and right rotating unstable modes of a straight vortex have the same growth rate, the resulting configuration 
is determined by the cross coupling coefficient between these modes, which is a function of the parameters b, c. 

The symmetry between left and right rotating helices can be broken by applying an additional twist to the straight 
filament. Studies of twisted filaments were performed by Rousseau et al (1998) and by Nam et al (1998). The 
simulations revealed stable helices and a secondary super-coiling instability. 

Nam et al (1998) performed linear stability analysis of a straight filament with twist. It was shown that the twist 
reduces the domain of stability for the straight filament. 

The question arising in this context is how to prepare a vortex with twist. One way to proceed is to create an 
inhomogeneity close to the axis of the vortex filament. The inhomogeneity will locally change the frequency of the 
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vortex and will result in persistent twist. This situation was realized experimentally in a reaction-diffusion system by 
Mironov et al (1996). 



VI. GENERALIZATIONS OF THE CGLE 

The CGLe is a minimal, universal model that cannot be further simplified. However, there are many ways to 
generalize it in order to include qualitatively new features (we are not concerned with additional terms that merely 
give quantitative corrections). The main trends of generalization can be associated with the different terms the CGLe: 

• generalization of the nonlinearity 

• generalization of the differential operator 

• generalization of the symmetry group 



A. Subcritlcal CGLe 

1. Small Amplitude Solutions in the weakly nonlinear case 

In this subsection we consider the CGLe with a destabilizing nonlinearity 

dtA = A + {l + ib)AA-{-l + ic)\A\'^A (80) 

The CGLe was first derived in this form for a physical system by Stewartson and Stuart (1971) in the context of plane 
Poiseuille flow where one has a (strongly) subcritical bifurcation. 

Spatially-homogeneous solutions of Eq. (80) with A{t = 0) = Aq diverge at finite time according the following 
expression: 



' ' [l + |yloP(l-e2«)]i/2 

Also, for c ^ the rate of phase winding diverges. When \A\ is a function of x, however, its behavior is more subtle. 

At first sight one expects blowup of the solutions which can be avoided by adding a fifth-order stabilizing term. 
It was suggested by Hocking and Stewartson (1972) and Hocking et al (1972), that for generic initial conditions the 
blowup does not occur in a considerable region of the parameter space (6, c). By considering the evolution of pulse-like 
solutions, the region in (6, c) plane was found where the solution remains bounded. Since |c|, which is the ratio of 
nonlinear dispersion over nonlinear growth, has to be sufficiently large, this may be called the "weakly subcritical 
case" . 

Then these results were apparently forgotten, and addressed again (independently) by Bretherton and Spiegel 
(1983) (in the limit c oo), and Schopf and Kramer (1991). They reproduced much of the results of Hocking and 
Stewartson (1972) and found stable periodic solutions of Eq. (80). The analytic work was supported by detailed 
simulations. The work was continued by Powell and Jakobsen (1993), Kaplan et al (1994a, b), Kramer et al (1995), 
and Popp et al (1998). Weakly subcritical Hopf bifurcations are found in convection in binary fluids (see Moses et al 
(1987), Heinrichs et al (1987), Kolodner et al, (1988, 1995, 1999)), and presumably in nonlinear optics (Powell and 
Jakobsen, 1993, Kramer et al, 1995). 

Hocking and Stewartson (1972), Bretherton and Spiegel (1983), and Schopf and Kramer (1991) have foimd bursts 
of two types depending on the relative signs of b, c (by analogy with the NLSe the case of 6c > will be called the 
focusing case and otherwise defocusing). 

Kaplan et al (1994a, b), Kramer et al (1995) proposed a simple physical mechanism, called the phase gradient 
mechanism, providing arrest of the collapse in Eq. (80) if |c| is sufficiently large and |6| not too large. The phase 
gradient effect manifests itself as a fast differential phase rotation resulting from the explosive burst amplitude increase. 
To understand this effect it is convenient to rewrite Eq. (80) in the variables A = \/i?exp[i6'] (see Eqs. (10)). For the 
sake of simplicity one considers the limit 6^1. Introducing the local wavenumber k = d^O Eqs. (10) reduce to 



dtR = (1 + i?2)i? + d^R - k^R 

i?2 



d,k = cd,R-' + dJ^^^^) (82) 
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From the second equation one sees that a gradient in R drives growth of |fc|, which in turn saturates R via the last term 
in the first equation. If |c| is sufficiently large, this effect overcomes the explosive growth manifest in the first term on 
the right-hand side of the first equation. Specifically, consider an initial pulse having a small, broad, constant plateau 
and decaying away at the edges Let fc = initially. At the linear stage of the instability and then in the forthcoming 
blow-up regime the amplitude will remain approximately constant inside the plateau region. At the boundaries sharp 
gradients of R will be formed. These gradients will act as sources for the generation of the phase gradient, i.e. k, in 
narrow regions. There the large value of k will saturate the blow-up. The net result will be two counter-propagating 
fronts representing the moving plateau boundaries. The front propagation speed will not be constant but it will grow 
at the blow-up stage. Thus, if |c| is sufficiently large, the pulse will be "eaten" by these fronts moving from the edges 
to the center. 

2. Strongly subcritical case 

The solutions c;onsidered above bifurcate from the trivial state supercritically (in the range of b, c where they remain 
bounded), in spite of the fact that the sign of the real part of the nonlinear term signalizes a subcritical bifurcation. 
Thus in that parameter range, but below threshold, the trivial state is the global attractor. The scenario is not 
changed qualitatively by the addition of stabilizing highcir-order terms. Outside this parameter range one needs (at 
least) quintic terms to saturate the explosive instability provided by the cubic term. The equation can be written in 
the form 

dtA = eA+{l + i6)A^ - (-1 + ic)\A\^A - (1 + id)\A\'^A (83) 

The finite-amplitude solutions persist stably with respect to amplitude fluctuations below threshold e < in a 
certain parameter range, where they coexist with the linearly stable trivial solution. There exist moving fronts and 
- surprisingly - stable localized pulses over a finite interval of parameters (Thual and Fauve, 1988). This clearly is 
a result of the nonvariational nature of Eq. (83). A study of existence, stability and selection of various solutions 
is given by van Saarloos and Hohenberg (1992). Localized perturbations around the trivial state can either decay 
(small e), or evolve into pulses (intermediate e < 0), or develop into fronts that invade a plane- wave state. In certain 
cases their velocity is selected by a "nonlinear" marginal stability criterion below some positive value of e and by 
linear marginal stability, as in the supercritical cubic CGLe, for larger e. In a parameter range with sufficiently small 
nonlinear dispersion there exists a class of fronts that can be expressed in the form (polynomial fronts) 

k = kN + eQ{R^ - R%), R' = eiR{R'^ - R%) (84) 

Very recently steady fronts that cannot be expressed in this form were found (CouUet and Kramer, 2001). They exist 
in particular in a parameter range where there are no polynomial fronts. Whereas the polynomial fronts are sources 
in their rest frame the new fronts are sinks. They move in a direction to generate the trivial state even for positive 
values of e, which can be understood from the phase gradient mechanism discussed above. They play a significant 
role in the dynamics of spatio-temporal chaotic states and for the formation of localized structures. 

Deissler and Brand (1994, 1995) and Akhmediev at at (2001) have studied periodic, quasiperiodic and chaotic pulses 
of Eq. (83) and its generalization. The interaction of these solutions in the framework of two coupled equations Eq. 
(83) shows that the result depends sensitively on the initial conditions. Convective and absolute instabilities in the 
subcritical CGLe are investigated by Colet al (1999). 

Deissler and Brand (1991) found in 2D localized particle-like solutions of Eq. (83). Moreover, one expects localized 
solutions possessing a topological charge coexisting with extended (conventional) spirals known for CGLe (Malomed 
and Rudenko, 1988). Recently, the properties of spiral waves and other localized solutions in the cubic-quintic CGLe 
were studied by Crasovan et al, 2001. 

B. Complex Swift-Hohenberg equation 

The complex Swift-Hohenberg equation in the form 

dtA = rA-{\ + ic)\A\'^A + ioAA - (fi A)^^ (85) 

where r <C 1 is the control parameter, a characterizes the diffraction properties of the active media, was derived 
asymptotically in the context of large aperture lasers (class A and B) with small detuning Q. between the atomic and 
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cavity frequencies (see Staliunas, 1993; Lega, Moloney and Newell, 1994, 1995). This equation is believed to be relevant 
also for oscillatory convection in binary fluids, however, it can not be derived asymptotically from appropriate Navier- 
Stokes equations. Eq. (85) is a generic equation in the vicinity of a codimension-2 bifurcation where the coefficient 
in front of the diffusive term is allowed to change sign (see. e.g. CouUet and Repaux, 1987). 

Clearly one should distinguish between the real Swift-Hohenberg equation (see, e.g. Cross and Hohenberg, 1993) 
and Eq. (85). The Swift-Hohenberg equation is a phenomenological model and can not be derived from the original 
equations . In contrast, the complex Swift-Hohenberg equation in the form (85) is derived asymptotically rigorous 
in the limit of ^ 0, i.e. in the so-called "long-wavelength" limit. In this limit the differential nonlinearities have 
formally higher order and, therefore, can be dropped. In this case the "wavenumber selecting term" {fl + A)'^A is 
just a small correction to the "diffraction term" iaAA. Thus, the complex Swift-Hohenberg equation can be treated 
as a perturbed CGLe. In contrast to the real Swift-Hohenberg equation, Eq. (85) has two independent wavenumber 
selection mechanisms: the first is related to the {fl + A)^ term, providing maximal amplification for the plane wave 
ex])[i{ujt + kx)] with the optimal wavenumber k = \/f2; the second is related to the selection of wavenumber by 
topological defects (e.g. spirals, holes) in the CGLe and relies on the diffraction iaAA. For this case the last term in 
Eq. (85) changes slightly the wavenumber selected by the defects (Aranson, Hochheiscr and Moloney, 1997). 

The stability of plane waves in Eq. (85) was studied in details by Lega et al (1994). A new feature is the zig-zag 
(transversal) instability of plane waves for wavenumbers away from the band center. 

In 2D Eq. (85) possesses, as the CGLe, topological defects in the form of spiral waves. Aranson, Hochheiser 
and Moloney, 1997 show that in the limit of small r these spiral waves undergo a core instability, leading to stable 
meandering. Another feature of Eq. (85) in 2D is the existence of domain boundaries between traveling waves with 
different orientation, usually called zipper states. Furthermore, in 2D, the question of wavenumber selection can be 
transformed into one of wavevector selection, since the domain wall can adjust the direction of ingoing or outgoing 
waves. The domain wall itself may no longer be stationary, but may move in a certain direction, if there is no reflection 
symmetry of the wave pattern with respect to the domain wall axis. The second spatial dimension (along the domain 
wall) opens the possibility for additional instabilities of the wall, as it was observed experimentally in convection in 
binary mixture by Moses et al (1987), La Porta and Surko (1997). Aranson and Tsimring (1995) have shown that near 
threshold the active (emitting waves) zipper states are alway unstable with respect to transversal undulation. The 
nonlinear stage of this instability leads to creation of a chain of topological point defects (spirals) which themselves are 
unstable. The latter appears to be analogous to the famous Kelvin-Helmholtz instability of a tangential discontinuity 
of shear flows. Passive (absorbing) zipper states turn out to be stable. 

C. CGLe with broken gauge invariance 

It is interesting to break the global gauge invariance of the CGLe. This corresponds in particular to a situation where 
a system undergoing a Hopf bifurcation with the frequency u)c is parametrically forced (modulated) at a frequency 
near uic/2. This leads to 

dtA = A{€ + iuj) + (1 + ih)AA - (1 + ic)\A\'^A + -yA* (86) 

Obviously, instead of global gauge invariance A Ae** one is left with the discrete symmetry A —A. Here (j is 
the detuning, 7 > is the amplitude of forcing and e describes the distance from the threshold of instability. 

1. From oscillations to bistability (e > 0) 

For e > 0; "'^^ depending on the values of the other parameters, Eq. (86) describes an oscillatory or bistable situation. 
In the latter case the system is in general (i.e. for b,c,u 0) of excitable nature. In the context of ferromagnets 
in a static magnetic fleld it describes domain walls separating two stable states (domains with opposite spins). The 
domain walls exhibit a transition involving a spontaneous breaking of chirality (in ferromagnets Ising walls become 



The Swift-Hohenberg equation is obtained by keeping only one nonlinear term. However, nonlinear terms involving also 
derivatives have formally the same order, and, strictly speaking, cannot be neglected 
^''Note that for e > 0, by rescaling t, r, u), 7, A in Eq. (86), e can be replaced by 1. 
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Bloch walls as the strength of crystal anisotropy is reduced, Lajzerowicz and Niez, 1978, 1979). In such an equilibrium 
situation the imaginary coefScients in (86) vanish and Eq. (86) can be cast into variational form: 

O.A=-^^ (87) 

where the "free energy" functional is of the form 

-e|A|2 + ^ + |VA|2 - l{{A*f+A^)\ dxdy (88) 



Equation (87) possesses stationary kink-like solutions connecting stable homogeneous equilibria A = ±-^1 + 7 for 
e = 1: 

Ai ^ ±^1 + 7tanh(Vl + 73;) (89) 

Ab = ± v'rT7tanh(v^a;) ± i . (90) 

cosh(v27a;) 

The first solution, called by Coullet et al (1990) Ising wall, is stable when 7 > 7c = 1/3 and the second solution, 
which is stable for 7 < 7^ and breaks chirality, is called Bloch wall. The order parameter A vanishes at the core of the 
Ising wall but not at the core of the Bloch wall. For 7 = 1/3 an exchange of stability (pitchfork bifurcation) occurs 
between these two solutions. 

Coullet et al (1990) investigated the behavior of Ising and Bloch walls in nonequilibrium conditions, when at least 
one of the coefficients b,c or lo is nonzero. They found that nonpotential effects in general lead to motion of the 
Bloch walls, and not of the Ising walls ( see also Sakaguchi, 1992; Coullet and Emilson , 1992; Mizuguchi and Sasa, 
1993; Chate et al, 1999). Coullet et al (1991) applied this concept to the description of the Ising-Bloch transition in 
ferromagnets in a rotating magnetic field. 

Frisch et al (1994) have studied Eq. (86) in 2D in the context of a homeotropically aligned nematic liquid crystal 
in a rotating magnetic field in the vicinity of the electric Fredericks transition. They have found that Bloch walls 
containing a defect separating the two variants assume the form of rotating spiral waves. Those spiral waves had been 
studied experimentally also by Migler and Meyer (1994). Frisch et al (1994) have shown that spiral waves in Eq. (86) 
combines properties of spirals in oscillatory (as CGLe) and excitable spirals (as in reaction-diffusion systems, Tyson 
and Keener, 1988). The problem of the spiral's frequency selection was solved by Aranson (1995) for small lo and 
b = c = . In particular, it was shown that the frequency is selected by the local curvature of the moving Bloch wall. 
For 7, w ^ 1 the spiral wave solution can be described in the phase approximation by the overdamped sine-Gordon 
equation (0 — arg A) 

dt(f) = (jo-') sin(2(/)) -t- A(/) (91) 

Korzinov et al (1992) have studied Eq. (86) in the context of periodically forced convection. Hanusse and Gomez- 
Gesteira (1994) considered it in the context of chemical systems. 
Frisch and Gilli (1995) as well as Coullet and Plaza (1994) considered the more general equation: 

dtA = A{1 + iuj) + (1 + i6)AA - (1 + ic)\A\'^A + 7A* + 70 (92) 

The term 70 is responsible for the effect of a tilted magnetic field. In contrast to Eq. (86) the spiral waves in the frame- 
work of Eq. (92) exhibit a diverse variety of behaviors from rigid rotation to meandering and even hypermeandering. 
The Ginzburg-Landau equation with more complicated forcing terms was studied by Gilli and Gil (1994). 



2. Parametric excitation of waves in the GLe 



Eq. (86) for e < is a phenomenological model of parametric excitation of surface waves in fluids (GoUub and 
Langer, 1999). In the absence of parametric driving 7A* the system always relaxes towards the trivial state A = Q. 
However, if the parametric driving exceeds the critical value 7c = (w — tb)/^/l + 6^, the trivial state ^1 = becomes 
unstable with respect to standing waves with wavevector = {hu -|- e)/(H- 6^) of arbitrary orientation. In this sense 
Eq. (86) is reminiscent of the Swift-Hohenberg equation (Cross and Hohenberg, 1993). 
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In 2D a variety of nontrivial static and dynamic states are possible. Theoretical studies of spatio-temporal chaos 
in Eq. (86) in the context of surface waves in fluids were conducted by Zhang and Vifials (1995). Spiral waves were 
studied experimentally and theoretically by Kiyashko et al (1996). 

Tsimring and Aranson (1997), Aranson et al. (1999) used Eq. (86) coupled to an additional field to describe pattern 
formation in a thin layer of vibrated granular materials in connection with experimental studies (Melo et al 1994, 
1995, Umbanhowar et al 1996). In this case the parameters 7 and u) can be associated with the amplitude and the 
frequency of external periodic driving. Depending on the values of 7 and lo, a variety of stable solutions ranging from 
localized oscillons and interfaces separating domains of opposite polarity to periodic stripes, squares and hexagons 
were found. The topology of the transition lines between different types of the solutions turns out to be in agreement 
with the experiments. 

D. Anisotropic CGLe in 2D 

In many physically relevant situations the 2D CGLe is essentially anisotropic (see Sec. I, class iii), i.e. it is of the 
form 

dtA = A+{l + ih)d^^A + (1 + ib2)dlA - (1 + ic)\A\^A (93) 

with bi 7^ &2- The equation in this from was studied by Weber et al (1991), Brown et al (1993), Roberts et al (1996). 
The stability of plane waves depends on the orientation of the wavevector. In particular, one has the situation when 
the wave is stable only in one direction. Under such condition Weber et al. (1991) found stable lattices of defects. 

New features of phase and defect chaos in the 2D anisotropic CGLe were found by Faller and Kramer (1998) and 
(1999). The phase-chaotic states exist in a broader parameter range than in the isotropic case, often even broader 
than in one dimension. They may represent the global attractor of the system. There exist two variants of phase 
chaos: a quasi-one dimensional and a two-dimensional solution. The transition to defect chaos is of intermittent type. 

E. Coupled Ginzburg-Landau Equations 

When both, the critical wavenumber Qc, and the critical frequency u>c are non-zero at the bifurcation (class iii) with 
reflection symmetry, Sec. I), the primary modes are traveling waves which in ID or in the presence of anisotropy are 
described by two coupled complex Ginzburg-Landau equations. The physical fields in the weakly nonlinear regime 

are of the form ^ Ab. cxp[—i{ujct — Qcx)] + ALCxp[—i{uct + Qcx)] + c.c, where Ar and Al are the complex amplitude 
of right/left traveling waves. In ID the coupled CGLes are given by 

dtAR + sO^Ar = Ar + {1 + ib)dlAR - ((1 + ic^\Ar\'' + (1 + «rf),g|Ai|2) Ar 

dtAL-sd^AL=AL + {l + ib)dlAL-{{l + ic)\AL\'' + {l+id)g\AR\'')AL (94) 

where ,s is the linear group velocity, {1 + id)g \s the complex coupling coupling coefficient between the two modes 
(Cross and Hohenberg, 1993). 

In addition to the usual GLe parameters b and c one here has three relevant parameters s, d and g. Elaborate 
surveys of the various regimes occurring in Eqs. (94) are given by van Hecke et al (1999), Riecke and Kramer (2000), 
see also, Amengual et al (1996,1997), Neufeld et al (1996). 

The case of strong suppression corresponds to g > 1. In this situation dual- wave solutions with Ar = Al are 
unstable. In contrast, single-wave solutions Ar ^ 0, Al = or vise versa can be stable. For g > 1 a. variety of source 
and sink solutions in counter-propagating waves were analyzed by Malomed (1994), Alvarez et al, (1997), van Hecke 
et al (1999). 

F. Complex defects in vector Ginzburg-Landau equation 

The vector CGLe (VCGLe) can be viewed as a particular case of Eq. (94) for s = 0: 

dtA± =A± + {1 + i6)VM± - (1 + ic) (|A±|2 + (1 + id)g\A^\^) A± (95) 

The problem of nonlinear dynamics of a complex vector field arises most naturally in the context of nonlinear optics, 
where the order parameter is the electric field envelope in the plane normal to the direction of propagation; the fields 
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A± can be identified with the two circularly polarized waves of opposite sense (Gil, 1993; Pisnien, 1994a, b, Haelterman 
and Sheppard, 1994, San Miguel, 1995, Hernandez-Garcia et al, 1999, Hoyuelos et al, 1999, Byryak et al, 1999). 

A distinguished feature of the VGLE is the possibility of a transition between two "phases" , which can be charac- 
terized by either "mixing" or "separation" of two "superfluids" . Defects (vortices) can exist in both "supcrfluids" , and 
transitions between alternative core structures are possible (Pismen, 1994, 1999); in this sense the real case (VGLe) 
could be viewed as a toy model for the 9-component description of superfluid ^He, dressed down to two components. 
Pismen (1994)a,b introduced the notation of "vector" and "scalar" defects in the VGLe. Vector defects have topo- 
logical charge (and, therefore, zeros of A± ) in both fields, whereas scalar defects have non-zero charge only in one of 
the fields A±. 

Simulations of the VCGLe by Hernandez-Garcia et al (1999,2000) have shown spiral wave patterns with an ex- 
ceptionally rich structure where both separated (but closely packed) "scalar" defects in the two fields and "vector" 
defects with a common core could be seen. 

A particularly intriguing possibility, suggested by Pismen (1994)a, is the formation of a boimd pair of defects in 
the two fields, i.e. a vortex "molecule" with dipole structure. Aranson and Pismen (2000), have shown that such 
a "molecule" requires complex coefficients. Analytical calculations were conducted in the limit of small coupling g 
between two complex fields. As was shown, the interaction between a well-separated pair of defects in two different 
fields is always long-range (power- like), in contrast to the interaction between defects in the same field which falls off 
exponentially as in a single CGLE (Aranson, Kramer, Weber, 1993a). In a certain region of parameters of the VCGLe 
stable rotating bound states of two defects - a "vortex molecule" - are found. 

G. Complex oscillatory media 

Two-dimensional oscillatory media exhibit a wide variety of wave phenomena, including spiral waves, phase and 
defect turbulence etc. The CGLe describes the dynamics of oscillatory media in the vicinity of a primary Hopf 
bifurcation. Brunnet et al (1994), Goryachev and Kapral (1996) have shown that spiral waves may also exist in 
systems with more complex local dynamics, e.g. period-doubling bifurcations or chaos. In this situation the rotational 
symmetry of spiral waves may be broken by line synchronization defects. 

Goryachev et al (1999) studied transitions to line defect turbulence in complex oscillatory media supporting spiral 
waves. Several types of line defect turbulence were found in a system where the local dynamics is described by a 
chaotic Rossler oscillator. Such complex periodic spirals and line-defect turbulence were observed experimentally in 
chemical systems by Park and Lee (1999). 

VII. CONCLUDING REMARKS 

In this work we attempted to overview a wide variety of dynamic phenomena described by the CGLe in one, 
two and three dimensions. The CGLe exhibits in many respects similar behavior in all dimensions, e.g. active and 
passive defects, distinct chaotic states, convective and absolute instabilities etc. Surprisingly, quantitative aspects of 
instabilities and transitions are different. In particular, the core (acceleration) instability of ID Nozaki-Bekki holes, 
2D spirals and 3D vortex filaments has different manifestations: Nozaki-Bekki holes undergo a stationary instability, 
spirals exhibit unsaturated Hopf bifurcations and 3D vortex filaments show a supercritical Hopf bifurcation. One 
observes a general trend in the region of occurrence: ID defects have the smallest stability domain, 2D spirals are 
stable in a much wider range of parameters and the stability domain again decreases for 3D vortex filaments. 

The unique combination of all these features in one equation stimulates continuous interest in this topic in a broad 
scientific community. The insights obtained from the CGLe over the last decades had an enormous impact on the 
physics of non-equilibrium systems, pattern formation, biophysics etc. and will be usefull for further progress in the 
physics of complex systems. 

Let us discuss briefly some open problems in the CGLe world. 

• Description of turbulent states in all dimensions. Although considerable work was done to identify the stability 
limits and transition lines between various turbulent states in the CGLe, surprisingly little is known on the 
statistical properties of these states. The main obstacle is the lack of appropriate analytical tools, since the 
traditional methods of statistical physics are not suitable for the description of spatio-temporal chaos. 

• The structure and statistical properties of the vortex glass. The "glassy" properties of this state (such as 
power-like decay of correlations, hierarchy of the relaxation times etc.) are not exposed yet. 
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• The revolutionary development in computers will make possible detailed investigation of the dynamics in 3D 
CGLe. 

Hopefully, the last two question will be elaborated on during the next decade. However, the problem of spatio- 
temporal chaos and turbulence will possibly require considerable time and effort. 
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FIG. 1. Absolute stability limits Qai,2 for four cuts in the b, c plane. Also included are the convective (=Eckhaus) stability 
limits that separate the stable (light shadowing) from the convectively unstable (dark shadowing) regions. 
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FIG. 2. Phase diagram of the ID CGLe. BFN = Benjamin- Feir-Newell instability line up to which one has convectively stable 
planc-wavc solutions; AI= absolute instability line up to which one has convectively unstable plano-wavo solutions; DC = 
boundary of existence of DC towards small |c|. The other lines pertain to standing Nozaki-Bekki hole solutions: CI=core 
instability line, EH = convective (Eckhaus) instability of the emitted plane waves, AH = absolute instability of the emitted 
plane waves, MOH = boundary between monotonic and oscillatory interaction. 
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FIG. 3. Acceleration instability in the perturbed (cubic) CGLe. Comparison of the reduced acceleration dtv/v from theory 
(full line, Eq.(43)) and simulations (squares) for b = 0.5, c = 2.0, |d| = 0.002 and varying phase arg{d), Stiller et al (1995)a,b. 
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FIG. 4. Simulations showing the velocity of the hole center v = v{t) of interacting hole-shock pairs (periodic boundary 
conditions). In Figures (a),(b),(c) the CGLe parameters were b = 0.5, c = 2.3, d' — +0.0025 (i.e. far away from the core 
instability line), (a) relaxation into a constantly moving solution for period P = 48.4. (b) selected final state with (almost) 
harmonic oscillating velocity for period P = 37.0. (c) selected final state with anharmonic oscillating velocity for period 
P — 40.0. (d) CGLe parameters b — 0.21, c = 1.3, d' = —0.005 (near the core instability line), d" = 0; state with oscillating 
velocity (including change of the direction) for period P = 50. 
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FIG. 5. Snapshot of the modulus |A| = |A(a;)| of a stable uniformly moving hole-shock pair in a simulation for b = 0.5, c = 
2.3, d' = +0.0025. The solution is space periodic with period P = 48.4. 
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FIG. 6. (a) Example of the bifurcation diagram of the quasiperiodic solutions for c — —2.0, p = 2tv/50 (see text). The inset 
illustrates the drift pitchfork bifurcation {v = branch not shown beyond bifurcation) . (b) QPS profiles at lower (full circle) 
and upper (open eirclc) branch, and at the saddle- node (SN) (triangle). HB and DN mean Hopf bifurcation and drift pitchfork 
bifurcation, Brusch et al, 2000 
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FIG. 7. Phase turbulence observed for b = —1 and c = 1.333. Part (a) is a spatio-temporal plot of A, part b shows snapshot 
Chate, 1994. 
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FIG. 9. Space-time plots (over a range of 512 x 1000) of |^| showing zigzagging holes near the transition to plane waves for 
b = —0.6, c = 1.2 (a); evolution of a homoclons in a background state with wavenumber 0.05 for b = — 0.6, c = 1.4, space x 
time = 512 x 250 (b), van Hecke, 1998. 
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FIG. 10. Space-time plot of the zeros of the local phase gradient dx<f> = (with A — \ A\ exp[i</)|]) for the perturbed cubic CGLe. 
The parameters are b = 0.9, c = 5.55 (i.e. slightly below the absolute instability line for the plane waves AH of the Figure 2) 
and d' = 0.005 (upper figure), d' = —0.005 (lower figure), d" = 0. Periodic boundary conditions and small amplitude noise 
as initial couflitions. Isolated lines can be identified as holes and shocks (alternating) which are spontaneously formed out of 
the chaotic state. In their neighborhood one has approximately plane waves. For decelerating d' the holes suppress the chaotic 
state after some transient (lower figure) and one is left with a "ID vortex glass" , while in the accelerating case (d' > 0, upper 
figure) , the chaotic state persist. 
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FIG. 11. Image of the amplitude |^| and contour lines of the phase 4> = axgA = 0, tt for single-charged spiral solution in a 
domain with periodic boundary conditions. Note edge vortex in the corner. Image is coded in the grey scale, maximum of the 
field corresponds to black, minimum to white. 
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FIG. 12. Stability limits of spiral wave solution. The boundary of convcctivc {EI) and absolute instability [AI) for the 
waves emitted by spiral are also shown (for explanation see (Aranson et al, 1992)), Bound states exist right of OR line 
(c — 6)/(l + be) = 0.845. BFN indicates Benjamin- Feir-Newell limit 1 + bc, L shows limit of 2D phase turbulence, and line T 
corresponds to transition to defect turbulence for random initial conditions (Chate and Manneville, 1996). 
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FIG. 13. Stability limits of spiral wave in the large b limit, e = 1/6. Here CI is the core stability limit e = Ec (unstable to the 
left), below EI is the Eckhaus unstable region, ST designates the transition line to strong turbulence (see text), below OR is 
the oscillatory range (Aranson et al, 1994). 
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FIG. 14. The distance L between two defects of opposite polarity that approach each other on a straight line is plotted for 
different Q versus time T (Bodenschatz et al, 1991, For comparison the experimental data of Braun and Steinberg (1991) are 
included. The different symbols denote different distances e from the threshold (circles: e = 0.02; squares e = 0.04; diamonds 
e = 0.005) 
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(a) b) 
FIG. 16. (a) Radial velocity (vx) and tangential velocity {vy) vs spiral separation 2X for 6 = 0, c = 1 for oppositely charged 
pair. Solid / dashed lines show radial and tangential velocities obtained from solution of Eq. (59) correspondingly. Symbols (o 
is radial and <> is tangential velocities) arc obtained from the simulations. Inset: Equilibrium distance 2Xe (solid line) given by 
Eq. (61) as function of c for 6 = 0, o and A correspond to numerical results for oppositely and like-charged pairs respectively, 
b) Tangential velocity Vy at the equilibrium distance for 6 = 0, • and A correspond to to numerical results for oppositely and 
like-charged pairs respectively. Solid and dashed lines show theoretical results obtained using nonlinear and linear treatments 
of the shock correspondingly. 
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(a) b) 
FIG. 17. (a) The tangential and radial velocities of the spiral versus the distance X from the inhomogeneity ^{r) = bo exp(— — ) 
for 6 = 0, c = 1 and bo = 0.3, a = 20. From the numerical simulations one finds the radius of the first stationary orbit ro »i 9.8; 
(b) The limiting orbits of the spiral core for different values of inhomogeneity "strength" 6o = 0.1, 0.2, 0.3 in the "large b limit" 
for £ = l/b = 0.2, c = 1, fT = 20 




FIG. 18. The dependence of radial velocity {vx) and the frequency increment C = d{dt'p)/d(p at = for 6 = 0, c = 1. 
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FIG. 19. Evolution of an unlike-charged spiral pair into a stable antisymmetric state, b = 0, c = 1, = Ly = 100, time lapse 
between snapshots = 50 




FIG. 20. Vortex glass in the convectively unstable range, 1024 x 1024 lattice, periodic boundary conditions (Chate and Man- 
neville, 1996). Large spirals nucleating in turbulent sea, c = 0.75, 6 = — 2 (left) and developed vortex glass for c = 0.7, b = —2 
(right). 
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FIG. 23. The number of defects versus time for four different values of c and e = 0.025. Parameters of the simulations are: the 
domain of integration 150 x 150, the number of Fourier harmonics 256 x 256, the time step 0.02. 




FIG. 24. Snapshots of \A{x, y)\ in the intermittency range for four consecutive times (c = —0.5, e = 0.025; black: \A{x, y)\ = 0, 
white: |j4(a;,j/)| = 1). The times corresponding to a,b,c,d are indicated in previous Figure 23. 
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FIG. 25. Dynamic of vortex rings in the NLSe. (a) In-plane collision of two rings at 90° incidence, seen form the side; (b) 
in-plane collision of two rings at 120° incidence, seen from above, Koplik and Levine (1996) 
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FIG. 26. Nucleation and expansion of vortices in the GLe with background k — 0.4. The supercurrent js is given by js = 
k{l — k^). Shown are the 3D isosurfaces of \A\ = 0.4. (a-c) No thermal fluctuations, images are taken at times t = 36,48,80. 
(d-f). Amplitude of thermal fluctuations 7/ = 0.002, t = 24,48,80. (Aranson, Kopnin, Vinokur, 1999) 



59 




FIG. 27. The growth rate ReX{k) as function of A; for e = 0.02, c = .1. Solid hne is the theoretical result for fe <§; 1, dashed line 
with symbols the result of numerical solution of 3D CGLe. Inset: blow-up of small k region. (Aranson and Bishop, 1997) 



FIG. 28. Instability of a straight vortex filament. 3D isosurfaces of = 0.1 for e = 0.02, c = —0.03, shown at four 

times: 50 (a), 150 (b), 250 (c), 500(d). Similar dynamics is observed also for larger value of e (Aranson and Bishop, 1997). 
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FIG. 29. The dependence of filament length L on time. Solid line corresponds to e = 0.02, c = —0.03; dashed line corresponds 
to e = 0.02,c= -0.5. 




FIG. 30. Two snapshots of 3D isosurfaces of \A\ taken in the regime of spatio-temporal intermittency, e = 0.02, c = —0.5. (a) 
Left image corresponds to t « 620 for Fig. 29; (b) right image corresponds to t w 740. (Aranson, Bishop, and Kramer, 1998). 
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FIG. 31. Sequence of snapshots demonstrating the evolution of a vortex ring for e = 0.2 and c = 0.2. 




FIG. 32. Stability limits in three-dimensions (solid line) and two dimensions (dashed line), obtained from linear stability 
analysis. Symbols represents the limit of the two-dimensional instability, obtained by direct numerical simulation of the CGLe 
(Aranson, Kramer, Weber, 1994). Vortex lines and two-dimensional spirals are stable to the right of the respective lines. 
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